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ABSTRACT 


A  two  dimensional  boundary  layer  model  is  developed  for  the  study 
of  diurnal  cycles  in  boundary  layer  convergence  over  synoptic  scale 
sloping  terrain  and  thermally  inhomogeneous  surfaces.  To  properly 
resolve  the  terrain  configuration, a  local  terrain  following  coordinate 
system  is  developed.  The  local  coordinate  system  has  characteristics 
which  are  superior  to  the  terrain  following  coordinate  systems  commonly 
used.  Truncation  errors  are  significantly  reduced  in  the  local  system 
so  variations  in  the  flow  field  due  to  physical  forcing  are  easily 
discernable  from  those  due  to  truncation  error.  Also,  it  appears  that 
the  local  system  can  resolve  coarser  terrain  configurations  under 
numerically  stable  conditions. 

The  model  exhibits  a  sensitivity  to  latitude.  The  phase  of  the 
maximum  vertical  velocities  forecast  by  the  model  occurs  six  hours 
later  at  22°N.  vs.  38°N.  The  boundary  layer  model  is  apparently  not 
very  sensitive  to  radiative  flux  divergence  calculations  but  did  show 
a  marked  reaction  to  soil  type.  Both  the  amplitude  and  phase  of  the 
maximum  vertical  velocities  are  modified  by  changes  in  soil  type. 
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CHAPTER  1 


INTRODUCTION 

The  diurnal  variability  of  precipitation  over  the  United  States  was 
first  documented  by  Kincer  (1916)  wf«o  used  summer  rainfall  data  from  175 
stations.  He  noted  that  at  many  stations  in  the  central  and  north  cen¬ 
tral  United  States  more  rain  fell  at  night  than  during  the  day.  Wallace 
(1975)  conducted  a  comp-'eitensue  study  of  the  amplitude  and  phase  of  the 
diurnal  cycle  of  precipitation  over  the  United  States.  He  also  noted 
that  over  the  central  part  of  the  United  States  in  summer  thunderstorms 
have  a  maximum  frequency  of  occurrence  around  midnight  with  most  other 
forms  of  precipitation  occurring  later  in  the  night.  This  is  interest¬ 
ing  since,  in  general,  thunderstorms  not  associated  with  frontal  acti¬ 
vity  have  a  tendency  for  maximum  occurrence  in  the  afternoon  hours  when 
lapse  rates  may  be  steep  due  to  surface  heating. 

Hewson  (1937)  proposed  that  radiative  cooling  from  cloud  tops  could 
explain  the  nocturnal  thunderstorm  maximum  over  the  central  United 
States  in  summer.  Means  (1944)  attributed  the  nocturnal  maximum  to  low 
level  warm  advection.  Although  both  processes  may  act  to  destabilize 
the  atmosphere  at  night,'  it  has  not  been  shown  that  these  processes  are 
strong  enough  to  counter  the  stabilizing  influence  at  night  of  sensible 
heat  flux  from  the  ground. 

It  has  been  recognized  that  over  much  of  the  central  portion  of 
the  United  States  during  summer  the  atmosphere  is  convectively  unstable 
24  hours  a  day.  Bleeker  and  Andre  (1951)  suggested  that  under  unstable 


conditions  even  a  small  diurnal  cycle  in  the  boundary  layer  convergence 
and  the  accompanying  rising  motion  at  the  top  of  the  boundary  layer  may 
be  sufficient  to  control  the  timing  of  convection. 

Many  investigators  have  related  the  duirnal  cycle  in  boundary  layer 
convergence  in  the  central  United  States  to  the  occurrence  of  a  low 
level  jet  in  the  southern  Great  Plains.  Bonner  (1968)  determined  the 
maximum  frequency  of  occurrence  of  the  southerly  low  level  jet  to  be  in 
western  Oklahoma  with  the  axis  of  maximum  frequency  oriented  northeast¬ 
ward  into  northeastern  Kansas  and  northwestern  Iowa.  A  statistical 
study  by  Pitchford  and  London  (1962)  revealed  a  correspondence  between 
the  mean  axis  of  the  low  level  jet  on  127  summer  days  and  the  line  of 
maximum  nocturnal  thunderstorm  occurrence. 

Bonner  et  al  (1968)  used  kinematically  computed  vertical  veloci¬ 
ties  for  10  southerly  low  level  jets  and  found  ascending  air  of  a  sy¬ 
noptic  scale  in  the  downstream  portion  of  the  jets.  Bonner  concluded 
that  the  downstream  portion  of  the  jet  should  be  a  preferred  region  for 
the  formation  of  nocturnal  thunderstorms.  Means  (1952),  Curtis  and 
Panofsky  (1958),  and  others  have  presented  arguments  linking  the  low 
level  jet  and  its  associated  vertical  velocity  field  with  the  occur¬ 
rence  of  convective  activity. 

Many  theories  have  been  advanced  as  to  the  cause  of  the  low  level 
jet  in  the  Great  Plains.  Blackadar  (1957)  attempted  to  show  that  the 
low  level  jet  arose  from  an  inertial  oscillation  of  the  ageostrophic 
wind  vector  as  the  air  near  the  top  of  the  friction  layer  was  decoupled 
from  the  air  below  by  the  formation  of  a  nocturnal  inversion.  Wexler 
(1961)  proposed  that  the  jet  formed  as  a  result  of  the  northward  de¬ 
flection  by  the  Rocky  Mountains  of  a  shallow  layer  of  air  flowing 
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westward  across  the  Gulf  of  Mexico.  Bonner  (1968)  found  that  the  diur¬ 
nal  oscillation  of  the  wind  in  the  boundary  layer  was  important  in 
producing  a  southerly  low  level  jet.  He  also  found  that  the  synoptic 
conditions  which  favored  large  scale  jet  formation  were  those  which 
contribute  to  a  strong  pressure  gradient  across  the  Great  Plains  with 
a  smooth,  uninterrupted  flow  of  air  northward  from  the  Gulf  of  Mexico. 
These  synoptic  conditions  are  common  over  the  southern  Great  Plains 
during  summer. 

Several  models  of  the  boundary  layer  have  been  developed  in  order 
to  gain  further  understanding  of  the  diurnal  oscillation  of  the  boun¬ 
dary  layer  wind  and  the  related  low  level  jet.  The  effect  of  varying 
eddy  viscosity  in  a  constant  geostropic  field  was  modeled  by  Buajitti 
and  Blackadar  (1957)  and  Ooyama  (1957).  Holton  (1967)  developed  a 
model  which  included  the  effect  of  a  diurnal  temperature  cycle  over 
sloping  terrain  with  constant  eddy  viscosity.  Bonner  and  Paegle  (1970) 
developed  a  model  which  included  both  effects.  None  of  these  models, 
however,  were  capable  of  predicting  supergeostrophic  flows  as  strong 
as  those  observed  by  Bonner  et  al  (1968). 

All  of  the  models  mentioned  above  assummed  horizontally  uniform 
flow  which  is  a  valid  assumption  for  flows  where  the  Rossby  number  is 
much  less  than  one.  Bonner  et  al  (1968)  found  that  for  the  10  souther¬ 
ly  low  level  jets  examined,  the  average  horizontal  wind  speed  gradient 
approached  f/2  to  the  right  of  the  jet  core  (Rg=l/2)  and  f  to  the  left 
of  the  jet  core  (RQ=1).  Therefore  it  is  apparent  that  a  numerical 
model  of  the  boundary  layer  should  allow  horizontally  non-uniform  flow 
in  order  to  better  predict  the  strong  supergeostrophic  flows  that  are 
observed. 


In  1972  Paegle  developed  a  numerical  model  of  the  boundary  layer 
which  incorporated  the  effects  of  horizontally  non-uniform  flow  and 
diurnal  oscillations  of  the  eddy  stresses  and  geostrophic  wind  field. 
Paegle  and  Rasch  (1973)  conducted  exhaustive  tests  of  this  model  for 
idealized  flows.  Their  solutions  above  the  terrain  slopes  of  the 
Great  Plains  differed  significantly  from  the  steady,  linear,  Ekman  so¬ 
lutions.  The  model  was  capable  of  producing  strong  low  level  southerly 
jets  in  the  western  portions  of  moderate  anticyclones  over  very  slight 
terrain  slopes.  The  areas  just  downstream  of  the  jets  were  character¬ 
ized  by  boundary  layer  convergence  and  rising  motion.  It  appeared  that 
the  greatest  single  contributing  factor  to  the  nocturnal  jet  and  boun¬ 
dary  layer  convergence  pattern  was  the  Great  Plains  topography  with  the 
eddy  stress  oscillation  of  lesser  importance. 

Paegle  and  McLawhorn  (1973)  used  the  model  to  investigate  the  re¬ 
lationship  between  nocturnal  thunderstorm  activity  and  boundary  layer 
convergence.  An  analysis  was  made  of  11  Great  Plains  nocturnal  thun¬ 
derstorm  occurrences  that  had  no  obvious  synoptic  scale  support.  The 
time  and  space  phasing  of  the  vertically  integrated  boundary  layer  con¬ 
vergence  forecast  by  the  model  was  compared  with  the  time  and  space 
phasing  of  the  11  thunderstorm  occurrences.  The  results  of  the  inves¬ 
tigation  were  very  promising  and  seemed  to  corroborate  the  hypothesis 
that  nocturnal  thunderstorms  might  often  have  their  birth  in  boundary 
layer  processes.  There  were  several  weaknesses  in  the  model,  however. 
The  formulation  of  eddy  stresses  was  rather  arbitrary,  second  horizon¬ 
tal  derivatives  of  the  flow  field  were  neglected,  and  there  was  an 
implicit  assumption  of  neutral  stratification. 

Paegle  (1978)  relaxed  certain  assumptions  inherent  in  the  model 


used  by  Paegle  and  Rasch  (1973)  and  Paegle  and  McLawhorn  (1973)  by  per¬ 
mitting  the  model  pressure  field  to  be  modified  by  the  convergence 
field.  The  resulting  system  contained  gravity  wave  modes  which  were 
excluded  in  the  earlier  model.  He  used  the  new  gravity  inertia  wave 
model  to  examine  slope  induced  diurnal  gravity  wave  convergence  fields 
over  the  United  States.  He  found  that  the  magnitude  of  these  fields 
may  be  sufficiently  large  to  completely  compensate  low  level  divergence 
fields  induced  by  free  tropospheric  effects.  Sensitivity  tests  conduc¬ 
ted  with  this  model  indicated  that  the  topographically  induced  boundary 
layer  convergence  fields  over  the  United  States  were  sensitive  to  dissi¬ 
pation,  to  boundary  layer  heating  and  its  vertical  extent,  as  well  as 
to  the  ambient  circulation. 

It  has  become  apparent  that  in  order  to  further  an  understanding 
of  the  physical  processes  and  simulation  requirements  of  those  boundary 
layer  motions  that  are  determined  by  the  duirnal  cycle  and  which  play 
an  important  role  in  convective  weather  a  particular  type  of  boundary 
layer  model  is  needed.  The  preliminary  studies  of  Paegle  and  Rasch 
(1973),  Paegle  and  McLawhorn  (1973),  and  Paegle  (1978)  have  suggested 
that  the  diurnal  cycle  of  boundary  layer  convergence  over  the  central 
United  States  is  highly  dependent  upon: 

(1)  The  previous  history  of  about  one  to  two  days  of  the  boundary 
layer  flow  and  low  level  forcings. 

(2)  The  ambient  large-scale  circulation. 

(3)  The  evolution  and  vertical  structure  of  the  boundary  layer 
thermal  wave. 


General  circulation  models  have  only  a  few,  if  any,  explicit  fore¬ 
cast  levels  in  the  boundary  layer;  therefore,  these  large  scale  models 


generally  lack  the  vertical  resolution  and  physics  to  account  for  items 
(1)  and  (3).  Many  new  mesocale  models  have  recently  been  developed, 
(e.g.  Anthes  et  al .  (1977),  Pielke  and  Mahrer  (1977),  Perkey  (1977) 
and  others).  These  models  have  more  explicit  resolution  of  the  boundary 
layer.  On  a  smaller  scale,  Deardorff  (1974),  Orlanski  et  al .  (1974), 
Welch  et  al.  (1977),  Yamada  and  Mel! or  (1974)  and  others  have  performed 
detailed  boundary  layer  forecasts  in  models  with  30  or  more  computation¬ 
al  levels  within  the  boundary  layer,  but  ignore,  or  minimize  troposphere 
interaction.  These  mesocale  models  can  usually  be  run  out  to  only  one 
day  or  less  because  of  expense  and  boundary  contamination.  Therefore 
it  is  not  likely  that  item  (1)  can  be  satisfied. 

It  is  likely  that  the  above  mentioned  models  do  simulate  boundary 
layer  convergence  fields  adequately  for  the  particular  phenomena  that 
have  been  studied.  However,  the  adequacy  of  their  approaches  remains 
to  be  demonstrated  for  the  diurnally  oscillating  case  over  synoptic 
scale  sloping  terrain  where  boundary  layer  convergence  fields  may  often 
dominate  large  scale  processes.  This  is  especially  true  for  the  noctur¬ 
nal  convection  over  the  central  United  States  that  is  characterized  by 
sharp,  strongly  supergeostrophic  nocturnal  jets  in  the  lowest  kilometer. 

The  intent  of  the  present  research  has  been  to  develop  a  boundary 
layer  model  with  sufficient  resolution  and  physics  to  define  the  inter¬ 
acting  thermal  and  momentum  boundary  layers  on  duirnal  time  scales 
above  synoptic  scale  sloping  terrain  and  thermally  inhomogeneous 
surfaces. 


CHAPTER  2 


EQUATIONS 


2.1  System  of  Equations 

The  formulation  of  the  governing  system  of  equations  is  based  upon 
the  following  approximations:  (1)  Boussinesq,  (2)  hydrostatic,  and 
(3)  two-dimensional  (no  variation  of  the  meteorological  parameters  or 
terrain  in  the  north-south  direction).  The  resulting  boundary  layer 
model  equations  in  a  terrain  following  coordinate  system  are: 

2.1.1  The  east-west  momentum  equation 
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2.1.2  The  north-south  momentum  equation 
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.av 


3V> 


2.1.3  The  moisture  equation 


',2) 


it 

2.1.4  The  thermodynamic  equation 

if.  +  uli  +  *21  =  q  +  J_  (K  ii\ 
at  ax  ^z  v  az  v  eaz' 

2.1.5  The  soil  heat  conduction  equation 


3Ts  a  3Ts 
_ 1  -  JL  (k  _ L1) 

at  az  v  s  az  1 


(3) 


(4) 


(5) 
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2.2  Meteorological  and  Mathematical  Symbols 

The  meteorological  and  mathematical  symbols  are  defined  as  follows 

2.2.1  O  represents  a  standard  atmospheric  value. 

2.2.2  (  )'  represents  a  deviation  from  the  standard  atmospheric 
value. 

2.2.3  u  is  the  west-east  velocity  component  of  the  wind  which  is 
parallel  to  the  west-east  slope  of  the  terrain. 

2.2.4  V  is  the  south-north  velocity  component  of  the  wind  which 
is  parallel  to  the  south-north  slope  of  the  terrain. 

2.2.5  W  is  the  vertical  velocity  component  of  the  wind  which  is 
perpendicular  to  the  X-Y  plane. 
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2.2.6  #  is  the  west  to  east  slope  of  the  terrain. 

2.2.7  K^,  KQ,  are  the  momentum,  thermal  and  soil  exhange 
coefficients. 

2.2.8  Q  represents  heating. 

2.2.9  P  and  p  are  atmospheric  pressure  and  density. 

2.2.10  T  and  are  atmospheric  and  soil  temperatures. 

2.2.11  Ty  and  e  are  virtual  and  potential  temperatures. 

2.2.12  R  and  Cp  are  gas  constants  for  dry  air. 

2.2.13  q  represents  specific  humidity. 

2.2.14  f  is  the  Coriolis  parameter. 

2.2.15  g  is  gravitational  acceleration. 

2.3  Application 

The  equations  may  be  applied  in  vertical -horizontal  coordinates 
in  a  non-orthogonal  sigma  coordinate  system  such  as  Gerrity  (1967). 
They  also  have  an  interesting  interpretation  in  an  orthogonal  terrain 
following  coordinate  system  (Appendix  A). 
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CHAPTER  3 


COORDINATE  SYSTEM 


3.1  Terrain  Following  Coordinate  System 

The  important  relationship  between  terrain  configuration  and  boun¬ 
dary  layer  convergence  was  established  in  Chapter  1.  It  is  apparent 
that  one  of  the  main  requirements  for  proper  simulation  of  boundary 
layer  convergence  is  a  terrain  following  coordinate  system  that  can  pro¬ 
perly  relate  the  terrain  configuration  to  the  governing  system  of  equa¬ 
tions.  To  satisfy  this  requirement  a  new  terrain  following  coordinate 
system  was  established.  A  complete  derivation  of  this  sytem  is  con¬ 
tained  in  Appendix  A. 

The  u  component  of  the  momentum  equation  written  in  the  local 
terrain  following  coordinates  is  listed  as  an  example. 


2 

cos  ® 


2r]ulf  *  t 


o 

1-sin  •sint<j> 


2  2 
tan  ®+sec  « 


-lyJLy.  _  cos41.  fv  = 
JwsT  cos® 


r _ 1  -I  1  3P.'  r  sin®  -1  £.'  +  1  -  JL  ( v  Ml 

"  1  c-in^cin^  ”  3X  rn«itanz*  +  l  «  9  C0S®  32  "32 

1-sin  ®sin  4*  p  cos  ®tan  4+1  p 


(11) 


where  $  is  the  south-north  slope  of  the  terrain.  The  terrain  slopes 
encountered  during  this  research  are  sufficiently  small  (the  maximum 
slope  is  about  .4853X10  radians)  to  simplify  the  coordinate  system. 
The  u  component  of  the  momentum  equation  written  in  the  simplified 
coordinates  is  : 


k 


i 


j*. 
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lr+  “w  +  *1?  -  fv  =  ■  -IT'  -  9  sin*  +  w  '"m!?' 


1  aP  * 


9U  \ 


(12) 


The  governing  equations  were  integrated  out  to  28  hours  using  two 
boundary  layer  models  which  are  identical  except  for  the  coordinate 
systems.  One  model  used  the  local  terrain  following  coordinates  while 
the  other  model  used  the  simplified  version.  There  was  no  perceptible 
difference  in  the  forecast  fields  at  any  time  during  the  28  hour  fore¬ 
cast  period.  Based  upon  the  results  of  this  test,  the  simplified 
terrain  following  coordinate  system  was  adopted  for  the  boundary  layer 
model.  The  governing  equations  of  Chapter  2  were  written  in  terms  of 
the  simplified  coordinate  system. 

3.2  Linearized  Pressure  Gradient  Force 

The  pressure  gradient  force  term  in  the  u  component  of  the  momen¬ 
tum  equation  is  linearized  to  separate  out  the  hydrostatic  vertical 
portion.  Therefore,  the  horizontal  pressure  gradient  force  is  no  long¬ 
er  formulated  as  a  small  balance  of  large  terms  over  sloping  terrain. 
This  avoids  major  numerical  truncation  errors  in  the  vicinity  of  moun¬ 
tains  that  are  common  for  most  types  of  terrain  following  coordinates 
such  as  the  frequently  used  sigma  coordinates. 

3.3  Transformation  of  Vertical  Coordinates 

Another  important  requirement  for  the  proper  modeling  of  boundary 
layer  convergence  is  the  correct  simulation  of  the  evolution  and  verti¬ 
cal  structure  of  the  low  level  thermal  field.  This  simulation  requires 
a  very  fine  vertical  grid  spacing  in  the  air  and  soil  near  the  air-soil 
interface.  To  provide  this  vertical  grid  spacing  while  retaining  some 


1 


12 

degree  of  economy  the  following  logarithmic  height  transformation  is 
used  for  the  atmosphere. 

Z  =  ZQ  ec/L  c  -  L  ln(Z/ZQ)  (13) 

where  Z  s  old  height  coordinate 

£  *  new  height  coordinate  j 

Zn  =  roughness  height  I 

u  • 

L  =  1  meter 

In  the  transformed  coordinate  (c).  logarithmic  variations  in  Z  become 
linear  in  £,  thereby  substantially  increasing  the  accuracy  of  the 
numerical  model  in  the  surface  layer.  The  u  component  of  the  momentum 
equation  written  in  the  transformed  height  coordinate  becomes: 


3u.  +  ..3U.  +  (Lx 
at  +  uiT+  Z' 


iu  .  f v  a  .  1  3P ' 

3£  ~  3X 

P 


-  ^  9  sin® 

p 


3£ 


(14) 


Analagous  forms  hold  for  equations  (2)  -  (4). 

A  guadratic  height  transformation  was  used  for  the  soil, 
o  -Z+Zr 


Z  =  -Ls^  +  Z0 


5S  *  (' 


-0J/2 
L„  ' 


(15) 


where  Z  =  old  depth  coordinate 


£$  *  new  depth  coordinate 
l%  =  .005  meters 

The  soil  heat  conduction  equation  written  in  the  transformed  height 
coordinate  becomes 

3Tn  -  ‘Ks  /  Ls  \3/2  3Ts  Ks  (  Ls 

,Lz+Z^  3?s  4L2  -Z+Z0  3, 


_Ss5(  Ls 
3t  4l2  -‘-■‘-q 


•) 


32T 
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CHAPTER  4 


NUMERICAL  ASPECTS 

4.1  Horizontal  and  Vertical  Grids 

The  horizontal  grid  for  the  boundary  layer  model  consists  of  31 
grid  points  which  extend  across  the  United  States  from  130°W  to  70°W 
at  38°N  latitude.  The  relationship  between  the  grid,  state  boundaries 
and  certain  geographical  features  is  depicted  in  Figures  1  and  2. 

Horizontal  grid  spacing  is  2°  longitude  or  175,246  meters.  Actual 
terrain  data  from  the  Air  Force  Aeronautical  Chart  and  Information  Cen¬ 
ter  (ACIC)  was  used.  This  data  was  properly  smoothed  by  multiple  appli¬ 
cations  of  a  single  nine-point  operator  as  described  by  Shuman  (1957) 
to  insure  that  all  variations  in  the  smoothed  terrain  could  be  resolved 
on  the  2°  horizontal  grid.  Terrain  slopes  were  computed  by  centered 
space  differencing. 

The  vertical  grid  consists  of  the  roughness  height  at  .01  meters 
plus  14  computational  levels  in  the  soil.  The  atmosphere  comprises  20 
computational  levels  in  the  transformed  grid  plus  15  equally  spaced 
levels  up  to  a  hieght  of  1520  meters.  Computational  levels  for  the 
boundary  layer  model  are  depicted  in  Table  1. 

4.2  Differencing  Schemes 

Implicit  differencing  with  respect  to  the  diffusive  terms  provides 
stable  solutions  with  time  steps  of  300  seconds.  as.75,  e=.25  were 
used  in  a  modified  form  of  the  Crank-Nicholson  scheme  as  discussed  in 


Relationship  between  the  horizontal  grid  and  state  boundaries. 
(1)  California,  (2)  Nevada,  (3)  Utah,  (4)  Colorado,  (5)  Kansas 
(6)  Missouri,  (7)  Illinois,  (8)  Kentucky,  (9)  West  Virginia, 


Figure  2.  Relationship  between  the  horizontal  grid  and  certain 

geographical  features.  (1)  Sierra  Nevada  Mts.»  (2)  Great 
Basin,  (3)  San  Juan  Mts.,  (4)  High  Plains,  (5)  Ozark  Plateau, 
(6)  Mississippi  River  Valley,  and  (7)  Appalachian  Highlands 
and  Mts. 
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Table  1.  Atmospheric  and  Soil  Computational  Levels  (Meters) 


Atmosphere 


Quadratic  Transformation  Logarithmic  Transformation 


-  Roughness  Height  - 


.0165 

.0272 

.0448 

.0739 

.1218 

.2009 

.3312 

.5460 

.9002 

1.4841 

2.4469 

4.0343 

6.6514 

10.9663 

18.0804 

29.8096 

49.1477 

81.0308 

133.5970 

220.2650 

Constant  aZ-86.668 

306.932 

393.599 

480.267 

566.934 

653.602 

740.269 

826.936 

913.604 

1000.270 

1086.938 

1173.606 

1260.274 

1346.942 

1433.610 

1520.278 
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Paegle  et  al.  (1976).  The  strong  physical  forcing  and  damping  in  the 
model  allows  use  of  centered  second  order  space  differencing  approxima¬ 
tions  for  the  advective  terms  without  incurring  non-linear  instabili¬ 
ties.  The  finite  difference  form  of  the  governing  system  of  equations 
written  in  the  simplified  terrain  following  coordinates  and  transformed 
height  coordinates  follows  with  A=?1.5  and  B= . 5 : 


4.2.1  The  east-west  momentum  equation 

n+1  n-1  n+1  ,  n+1  n+1  ,  n+1  9lin+l 

/Ui  ,K  '  ui,K>  ,  A,ui,K+l  '  ui,K-^  A  „  A,ui  ,K+1  +  ui  ,K-1  ' 

__  }  +  L  ^  __  j  +  2 


n-1 


n-1 


n-1  ,  n-1  9  n-1 

„  B,ui,K+l  "  ui ,K-1 \  n  B/Ui ,K+1  ui,K-l  -  *uik  x 
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-  p,r? 
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+  _J —  +  LLi  g  sine. 

pirK  pi,K 
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4.2.2  The  north-south  momentum  equation 

I  -  v 
TaT 


wn+1  .  v0*1  Vn+1  -  vn+1  vn+1  +  vn+1  -2vn+1 
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2At  }  L  V  Ti C  ;  U  21  2  '  ~ 
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r  ByVljjHl  n  B,Vi ,K+1  Vi,K-l  ^Vi  ,K> 

"  L  21  2a?  }  '  U  2*  a?2  ' 
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*  +  *  K,k 
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In  the  north-south  and  east-west  momentum  equations 
v  n  „  n 

r  _  /  L  \2,  m  i ,K+1  i ,K-1  ^  „  n  ,  LN 

c  -  (^}  ( - TZi - >  ■  S  i.k  y 
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j  n+1  T  n-1  T  n+1  T  n+1  T  n+1  ,T  n+1  ~T  n+1 

/  ‘  s  i  ,K~  s  i  ,K\  .  .  A,  s  i  ,K+l"  s  i,K-K  .  „  A, 's  i,K+l  s  i.K-T^'s  i,KN 

( - 2H - 5  +C  2*  2lc  }  +  D  2>  ^2  } 


T  n-1  T  n-1 
r  s  1  »K+1"  s  1  >K-1\  n 
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n-1  ,,  n-1 

i ,K-1 s  i  ,KN 
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(25) 


where 

-K. 


C  = 


and 


4L 


D  = 


i,K  ,  Ls  ,3/2  ,  / _ I  „Ks  i,K+l  ~  Ks  i,K-l 
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4.2.6  The  hydrostatic  equation 


Pi,K  =  Pi,K+l+(Pi,K+l'Rpi,KTvni,K)(^S_^)(  K  2L  1)Ac 

l  ,K 


4.2.7  The  continuity  equation 
Wi,K  =  Wi,K-l  '  ^  2L 


n  =  „n  A  /ZK  +  ZK-1 wUi+l ,K  ~  Ui-1,K> 

' (  2ax  ' 


(26) 


(27) 


(28) 


(29) 


In  the  above  finite  difference  equations,  the  subscripts  i  and  K 
denote  horizontal  and  vertical  grid  positions,  respectively.  The  su¬ 
perscript  n  denotes  the  time  step.  All  other  variables  and  mathemati¬ 
cal  symbols  were  defined  in  Chapter  2  with  the  exception  of  Z  which 
denotes  height  above  the  terrain  surface. 

The  coefficient  matrix  for  each  system  of  finite  difference  equa¬ 
tions  is  written  in  tridiagonal  from  and  solved  by  an  efficient  form 
of  Gauss  elimination  (Lindzen  and  Kuo,  1969).  The  complete  boundary 
layer  model  requires  approximately  3,100  seconds  of  Univac  1108  CPU 
time  for  a  60  hour  forecast  or  approximately  four  seconds  per  five 
minute  time  step. 
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4.3  Boundary  and  Interface  Conditions 

Values  for  forecast  variables  u,  vand  q  are  specified  at  the  rough¬ 
ness  height  (Z=Zq)  and  at  the  top  of  the  model  ( Z=1 520  meters). 

At  Z=Zq  u=v=w=0 

q=q0 

At  Z=1520  meters  u  and  v  are  specified 


P1  is  specified 

where  (  )q  =  initial  value.  Lateral  boundary  conditions  specify  a 
closed  domain:  ^  (u.v.q.e.P*)  «  0,  i.e.  total  reflection. 

The  temperature  field  required  special  treatment.  The  soil  and 
atmospheric  temperature  fields  are  linked  through  conditions  of  heat 
balance  and  continuity  of  temperature  at  the  roughness  height.  The 
surface  heat  balance  equation  is 

G  -  F  +  »CPK6  H  +  »5CSKS  ¥=lltI  =  Z  0  1301 

where  G  =  global  radiation 

F  =  net  longwave  radiative  flux 
ps,cs  =  so^  dens''ty  anc*  conductivity 
The  temperature  continuity  equation  is 

Tair  '  Ts  “  Z=Z0  <3’> 

In  finite  difference  form  the  surface  heat  balance  equation  becomes: 


1 

1 
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r-  .  H  r  1/  ^  f  1  *  1 

F  +  pi,0CPKe  i,Ol  Z,  -  zr 


L)  -  psCsKs(~ 


(’JjlLq" 

M013;ei,( 


s  1,-1 


)  (32) 


Because  of  the  low  conductivity  of  the  soil,  the  soil  temperature  at  a 
depth  of  approximately  one  meter  is  assummed  to  undergo  no  diurnal  vari¬ 
ation  over  the  forecast  periods  of  the  model ;  therefore,  the  lower 
boundary  condition  for  temperature  was  ts=tSq  at  z=--97  meters- 


4.4  Moisture  Constraint 

The  moisture  field  forecast  by  the  model  is  constrained  to  prevent 
supersaturation.  At  each  time  step,  the  newly  forecast  moisture  field 
is  converted  to  a  dew  point  field.  The  dew  point  field  is  compared  to 
the  newly  computed  temperature  field.  Any  dew  point  value  larger  than 
its  corresponding  temperature  value  is  set  equal  to  the  temperature 
value.  The  moisture  field  is  then  recovered  from  the  corrected  dew 
point  field.  The  expression  which  relates  the  moisture  field  (q)  and 


the  dew  point  field  (Tg)  follows: 


q  =  0.622  {  *  exp  [ 


597-0. 285(Td-273)  1  1 


<5Tf  -  r->]  )  (33) 


where  TD  =  dew  point.  A  derivation  of  the  above  expression  is 
contained  in  Appendix  B. 
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CHAPTER  5 

EXCHANGE  COEFFICIENT 


5-1  KEYPS  Formulation 

A  search  for  the  appropriate  formulation  of  exchange  coefficient 
included  the  investigation  of  a  turbulent  energy  closure  scheme  (Morin 
and  Yaglom,  1970),  the  Blackadar  formulation  for  a  neutral  atmosphere 
modified  for  non-neutral  conditions  (Yu,  1977),  and  the  KEYPS  formula¬ 
tion.  Numerical  instabilities  were  encountered  when  the  turbulent  energy 
closure  scheme  and  Blackadar1 s  formulation  were  structured  for  the 
boundary  layer  model  vertical  grid.  The  KEYPS  formulation  functioned 
well  in  the  model  environment  and  was  adopted  for  the  boundary  layer 
model . 

Assumptions  inherent  in  the  KEYPS  formulation  are  best  shown  by  a 
brief  derivation  of  the  KEYPS  expression.  The  turbulent  energy  equa¬ 
tion  is  written  in  the  following  form: 


db 

dt 


+ 


vv'b 


v  w 


3Z 


e 


(34) 


where  (~)  denotes  horizontal  averages.  The  terms  on  the  left  side  of 
the  equation  are  neglected  and  the  following  approximations  are  made: 


u  w 
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_KM  3Z 
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v  3V 

"*M  3Z 


w'0'  = 


»  36 
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(35) 


Multiplying  by  K^,  equation  (34)  becomes 

Kfr#2*(f!)2]  -  K"m  #  -  V  ■ 0 


(36) 
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Let  Y  =  MI  and  through  the  usual  dimensional  considerations  (e.g. 

Si9 

Plate,  1971)  let 

£'  W 


where  C  =  constant 

sl  -  length  scale 
Equation  (36)  now  becomes 

i/2 r  /  3U\2+/ 5v \2-i  £  Yj(2  30.  _  Si  -  q 

nL  iT  J  0  YSi  3z  (C£)4  u 


(37) 


Note  that 


U*  = 


-  A 


Si 


3Z 


surface  layer 


where 


|v )  =  S2  +  v2 

and  U*  is  friction  velocity.  Making  these  substitutions  equation  (37) 
becomes 


i,4  g  VZ  30 

U*  e  yKM  i?  • 


•4 


(c*) 


=  0 


(38) 


Although  strictly  valid  only  in  the  surface  layer,  equation  (38)  is 
used  for  the  entire  boundary  layer.  Solving  equation  (38)  for  Kj^ 
gives  our  final  result. 


Km  v  Hf)  *  (nc*)2  T  Hf]2  +  4U*>'/2}1/2 


(39) 
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Values  of  the  constants  are  y  =  14  (Lumley  and  Panofsky,  1964)  and 
C=1  which  insured  well  behaved  and  realistic  values  for  K^.  Selection 
of  C=1  also  insures  that  the  KEYPS  formulation  for  matches  similarity 
theory  under  neutral  stratification  in  the  surface  boundary  layer,  i.e. 
K^  «  U*kz  where  k  is  von  Karmon's  constant. 

The  KEYPS  formulation  projected  a  great  sensitivity  to  the  verti¬ 
cal  profile  of  the  potential  temperature  and  tended  to  produce  larger 
than  normal  values  of  in  the  upper  portion  of  the  boundary  layer. 

To  offset  the  latter  tendency,  was  allowed  to  attain  a  maximum  value 
below  a  specified  level  and  then  its  value  was  exponentially  decreased 
with  height  to  the  top  of  the  boundary  layer. 


5.2  Length  Scale 

The  length  scale  used  in  conjunction  with  the  KEYPS  formulation 
was  defined  by  Djolv  (1973)  as 


i  - 


W 


kZ 


*H<U> 


(40) 


P  7 

where  X=0. 00027  j  (Blackadar,  1962),  ^(p  is  the  Obukhov  profile  func¬ 


tion  for  momentum,  L  is  the  Monin-Obukhov  length,  and  G  is  the  geo- 
strophic  wind. 


5.3  Exchange  Coefficient  for  Heat 

The  KEYPS  formulation  provides  an  expression  for  the  exhange  co¬ 
efficient  of  momentum.  In  order  to  calculate  the  exchange  coefficient 
for  heat,  the  following  relationship  from  similarity  theory  is  used: 


(41) 


& 
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The  functional  forms  of  and  $  are  adopted  from  Carl  et  al .  (1973) 
and  Webb  (1970)  and  applied  as  in  Zdunkowski  et  al .  (1976).  The  sim¬ 
plifying  assumption  that  $  =<j>_  was  made.  The  KEYPS  formulation  assumes 

o  M 

that  the  ratio  of  K0  to  is  constant  with  height;  therefore  there  is 
the  implied  assumption  that  the  ratio  of  <t>  to  ij>M  is  constant  with 
height.  In  keeping  with  this  basic  assumption  of  the  KEYPS  formulation, 
the  ratio  of  <pQ  to  ^  was  calculated  in  the  surface  boundary  layer  and 
held  constant  with  height.  This  procedure  was  also  followed  by  Yu 
(1977)  with  favorable  results. 


CHAPTER  6 


RADIATION 


6.1  Radiative  Effects 

Radiation  appears  explicitly  in  the  model  as  the  heating  term  in 
the  thermodynamic  equation  and  as  global  and  terrestrial  radiation  in 
the  statement  of  heat  balance  at  the  roughness  height.  The  radiative 
effects  vary  throughout  the  boundary  layer  as  these  three  terms  act 
and  react  to  other  influences  simulated  by  the  model. 

6.2  Global  Radiation 

Global  radiation  is  calculated  as  in  Zdunkowski  et  al.  (1975).  A 
semi -empirical  formula  from  Phillips  (1962)  is  used  which  combines  the 
effects  of  diffuse  and  direct  solar  radiation.  The  formula  expresses 
the  global  radiation  in  terms  of  Linke's  turbidity  factor,  x^. 

G  =  J  (.755  -  .054  tL)  cosez  (42) 

where  J  =  solar  constant 

ez  =  solar  zenith  angle 

Equation  (42)  assumes  an  average  ground  reflectivity  of  20%.  The  Linke 
Turbidity  factor  is  defined  as 

tl  =  (eR  +  Sw)/8r  (43) 

where  br  =  complex  Rayleigh  scattering  coefficient 

3^  =  complex  water  vapor  absorption  coefficient 


J 


From  Feussner  and  Dubois  (1930) 


gR  =  -  1  n( . 907  nT018)  (44) 

where  m  =  air  mass. 

The  air  mass  is  normally  approximated  by  sec  6q.  A  more  precise 
expression  from  Paltridge  and  Platt  (1976)  was  used  in  the  model.  The 
expression  follows: 

«n  =  {[(^  cosez]2  +  2(JCjj£.)  +  1}1/2  -  (^)  cos0z  (45) 

where  r  =  earth's  radius 

p,P  =  Atmospheric  density  and  pressure  at  Zq. 
g  =  gravity 

McDonald  (1960)  provided  an  approximation  for  &w  where 

6W  =  (“■)  In  ( 1  ~ •  077 u * 30 )  (46) 

u  =  U1  *  m  (47) 

U-j  is  the  total  water  vapor  pathlength  of  the  atmosphere.  It  was  cal¬ 
culated  using 

00  p 

U,  =  /  qp  p-  dz  where  Pn  =  1013  mb.  (48) 

'  0  *0  u 

Infinity,  as  the  upper  limit  of  integration,  was  actually  the  top  of 
the  model;  therefore,  in  order  to  obtain  more  realistic  values  for  the 
water  vapor  pathlength,  the  top  of  the  model  was  extended  five  kilo¬ 
meters  with  a  constant  vertical  grid  separation  of  approximately  86 
meters.  The  resulting  radiational  model  is  6.52  KM  deep  with  94  com¬ 
putational  levels.  An  average  mid-latitude  sunnier  atmosphere 


i 


nil  r  ii>  m 
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(McClatchey  et  al . ,  1971)  was  Incorporated  in  the  upper  five  kilometers 
of  the  radiational  model  (see  Table  2). 


6.3  Flux  Divergence 

The  heating  term  in  the  thermodynamic  equation  is  expressed  in 
terms  of  the  flux  divergence. 


Q  = 


-1  ,JjKR/CP  3F 
PCp  P  3Z 


(49) 


where  F  =  infrared  flux 
PQ  =  1013  mb 

The  infrared  flux  divergence  (r—)  was  calculated  in  accordance  with 

6  Z 

Liou  (1979).  Flux  up  (F+(U))  and  flux  down  (F+(U))  are  calculated  at 
each  level  of  the  radiational  model.  F(U)  is  defined  as 


F(U)  *  F+(U)  -  F+(U)  (50) 

The  flux  divergence  is  the  difference  in  F(U)  between  layers  of  the 
model . 

The  flux  equations  for  the  radiational  model  are  written  in  terms 
of  U.  Since  the  model  atmosphere  is  considered  to  contain  only  water 
vapor,  U  is  the  optical  pathlength  of  water  vapor  and  is  calculated 
from  the  following  expression: 

U(z)  =  /  q(z)p(z) (~p^“)  dz  (51) 

0  *0 

where  Pg  *  1013  mb 

The  value  of  U  increases  from  zero  at  the  roughness  height  to  a  maximum 
value  at  the  top  of  the  radiational  model  as  depicted  in  Figure  3.  The 
flux  equations  are: 
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Table  2.  Mid-latitude  summer  climatology  data  used  in  radiational 
model . 


■ 

1 

eight(KM)  Pressure(mb)  Temperature (°K) 

■5  2 

Den:;ity(g/m  )  Water  vapor (g/m ) 

0 

1013 

294 

1.191E+03 

1.4E+01 

1 

902 

290 

1.080E+03 

9.3E+00 

2 

802 

285 

9.757E+02 

5.9E+00 

3 

710 

279 

8.846E+02 

3.3E+00 

4 

628 

273 

7.998E+02 

1.9E+00 

5 

554 

267 

7.211E+02 

1.0E+00 

6 

487 

261 

6.487E+02 

6.1E-01 

7 

426 

255 

5.830E+-2 

3.7E-01 

8 

372 

248 

5.225E+02 

2.1E-01 

9 

324 

242 

4.669E+02 

1.2E-01 

10 

281 

235 

4.159E+02 

6.4E-02 

Figure  3.  Vertical  variation  of  the  water  vapor  pathlength,  U,  and  its 
relation  to  Z. 
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F+(U)  =  oT4  ( U )  -  4oJ  tF(U-U' )  T3  ( U  * )  dTjK )  dU1 
0  u 

F+(U)  =  aT4(U)  -  oT4(U1  )tF(U1-U)  -  4a/  tF (U '  -U )T3 ( U ' )  dT-^[PdU-' 

U1 

(53) 

where  a  =  Stefan-Boltzman  constant 

tF  =  broad-band  flux  transmissivity 

=  U  at  the  top  of  the  radiational  model 
U'  =  U  at  an  intermediary  computational  level 
The  temperature  of  the  surface  and  the  air  immediately  overlying  the 
surface  are  assummed  to  be  the  same  in  expression  (52) 

The  broad-band  flux  transmissivity  is  calculated  in  terms  of  the 
broad-band  flux  emissivity  where 

tF  =  1  -  eF  (54) 

z  is  assummed  to  be  a  function  of  U  only  (there  is  a  slight  variation 
of  eF  with  temperature).  Broad-band  flux  emissivity  values  were  ob¬ 
tained  from  Stanley  and  Jarica  (1970).  These  values  include  the 
rotational  band,  window,  and  6.3  y  band  for  water  vapor.  Functions 
were  developed  to  express  emissivity  values  in  terms  of  u(u=log-|Q  U) 
for  use  in  the  model.  The  functions  are: 

US-3  eF  =  .0276u2  +  .282u  +  .7436 

-3<u*-l  eF  =  .0045u2  +  .1535u  +  .566 

-l<u*l  eF  =  .0515u2  +  .2395u  +  .605 

1 <u  eF  =  -.13u2  +  . 4823u  +  .5437 


(55) 


32 


Infrared  cooling  rates  calculated  by  the  radiational  model  were 
compared  with  cooling  rates  from  Roewe  and  Liou  (157C)  for  a  clear 
tropical  atmosphere  (McClatchey  et  al.,  1971).  The  results  are  depicted 
in  Figure  4.  The  low  cooling  rates  in  the  lower  troposphere  are  mainly 
due  to  the  failure  of  the  emissivity  values  to  fully  account  for  the 
water  dimer  effect  which  is  largely  responsible  for  water  vapor  absorp¬ 
tion  in  the  lower  troposphere  (Roewe  and  Liou,  1978). 

6.4  Solar  Zenith  Angle 

The  solar  zenith  angle  is  calculated  by  the  following  well  known 
expression  (Paltridge  and  Platt,  1976; z 

cosez  =  sin6  sine  +  coss  cose  costh  (56) 

where  ez  =  solar  zenith  angle 
6  =  declination 
e  =  latitude 
t^  =  hour  angle 

The  maximum  change  in  <5  in  24  hours  is  less  than  5°  so  a  single  value 
of  6  was  used  for  each  day.  The  expression  is  (Spencer,  1971): 

<5  =  .006918  -  .399912  cos9D  +  .070257  sineD 

-  .006758  cos29q  +  .000907  sin2@g 

-  .002697  cos39q  +  .001480  sin39p  (57) 

6q,  in  radians,  is  defined  in  terms  of  the  day  number,  dn,  which  ranges 

from  0  on  1  January  to  364  on  31  December. 

2ird 

eD  =  leT 


(58) 
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The  sun's  hour  angle,  t^,  is  measured  from  the  midday  position  and 
changes  360°  per  day  or  15°  per  hour.  The  radiational  model  contains 
two  internal  clocks  which  up-date  the  day  number  and  local  time  at  each 
grid  point. 

6.5  Numerical  Considerations 

The  number  of  mathematical  operations  required  for  the  radiational 
computations  was  at  least  an  order  of  magnitude  larger  that  the  number 
of  operations  required  for  the  remainder  of  the  boundary  layer  model; 
therefore,  the  radiational  computations  were  completed  every  23  time 
steps.  More  details  about  radiation  are  given  in  Chapter  7. 


CHAPTER  7 


SENSITIVITY  TESTS 


7.1  Coordinate  System 

7.1.1  Purpose 

The  purpose  of  this  test  is  to  give  some  indication  of  the  rela¬ 
tive  merits  of  the  local  terrain  following  coordinate  system  used  in 
the  boundary  layer  model  as  compared  to  a  standard  terrain  following 
coordinate  system  (as  in  Gerrity,  1967). 

7.1.2  Test  Description 

The  local  terrain  following  coordinate  system  and  the  coordinate 
system  used  by  Gerrity  (1967)  differ  only  in  the  manner  in  which  the 
pressure  gradient  force  is  expressed.  In  the  local  coordinate  system, 
the  pressure  gradient  force  is  written  in  terms  of  the  pressure  and 
density  deviations  from  standard;  whereas,  in  the  coordinate  system 
employed  by  Gerrity,  it  is  expressed  in  the  normal  manner.  Since 
pressure  and  density  deviations  from  standard  may  be  of  concern,  the 
coordinate  systems  were  tested  under  various  degrees  of  stratification. 
Care  was  taken  for  two  of  the  cases  to  initialize  the  model  in  such  a 
way  that  little  or  no  flow  would  be  generated;  therefore,  the  amount  of 
flow  that  was  generated  would  be  a  measure  of  the  truncation  error  in¬ 
herent  in  the  numerical  computation  of  horizontal  pressure  derivatives 
over  sloping  terrain  in  terrain  following  coordinate  systems. 

The  tests  were  catagorized  by  initial  atmosphere  as  follows:  (1) 
Condition  A,  (2)  Condition  B,  and  (3)  Condition  C.  Conditions  A  and  B 
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had  the  same  initial  atmosphere  except  for  the  type  and  strength  of 
temperature  inversion  (see  Figures  5  and  6).  The  model  atmosphere  was 
considered  dry  and  in  hydrostatic  balance  with  the  U,  V,  and  W  compon¬ 
ents  of  the  wind  field  equal  to  zero.  Global  radiation  and  radiative 
flux  divergence  were  omitted.  Friction  was  considered  constant  with 
Kj^  and  Kg  equal  to  one  meter  square  per  second.  The  model  atmosphere 
for  Condition  C  was  more  realistic  with  moisture,  global  radiation, 
and  friction  retained.  A  southerly  geostrophic  wind  of  10  m/sec 
existed  with  the  U  component  of  the  wind  field  still  equal  to  zero. 

The  V  component  of  the  initial  wind  field  increased  logarithmically 
from  zero  at  the  roughness  height  to  10  m/sec  at  50m  and  was  then  con¬ 
stant  with  height  to  the  top  of  the  model. 

The  boundary  layer  model  was  executed  using  each  coordinate  sys¬ 
tem  with  the  three  initial  atmospheres  for  a  48  hour  forecast  period. 

7.1.3  Results 

A  detailed  analysis  of  six  grid  points,  shown  in  Figure  7,  is  pre¬ 
sented.  The  grid  points  were  carefully  selected  to  insure  that  the 
absolute  value  of  the  terrain  slopes  associated  with  these  grid  points 
were  representative  of  the  range  of  slope  values  found  on  the  grid 
(Table  3).  The  generation  of  the  flow  field  during  each  test  is  repre¬ 
sented  by  graphs  of  the  maximum  U  component  vs  time  as  well  as  graphs 
of  the  integrated  vertical  velocity  at  the  top  of  the  model  vs  time 
for  each  of  the  selected  grid  points.  These  graphs  are  depicted  in 
Figures  8  through  35.  Averages  compiled  from  these  figures  are  con¬ 
tained  in  Tables  4  and  5. 

Data  contained  in  the  figures  and  tables  reveals  an  apparent 
superiority  of  the  local  terrain  following  coordinate  system  over  the 
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Figure  5 


CONDITION  A 


The  temperature  inversion  portion  of  Condition  A  used  to 
test  the  local  and  standard  terrain  following  coordinate 
systems. 


CONDITION  B 


The  temperature  inversion  portion  of  Condition  B  used  to 
test  the  local  and  standard  terrain  following  coordinate 
systems. 


coordinate  systems 


40 


Table  3.  Terrain  slopes  for  grid  points  used  in  the  coordinate 
system  tests. 


Grid  Point 

|  Slope  | 

5 

.004853 

14 

.004508 

15 

.002322 

27 

.001746 

25 

.000993 

18 

.000485 

standard  system.  In  most  cases  examined,  the  amount  of  flow  generated 
(|UmaXl)  in  the  standard  terrain  following  coordinate  system  increased 
as  the  terrain  slope  increased;  whereas  in  the  local  coordinate  system, 
the  flow  generally  increased  from  Condition  A  to  Condition  B.  These 
statements  also  apply  to  the  generation  of  the  vertical  velocity  fields 
in  the  standard  coordinate  system  but  only  for  the  larger  terrain 
slopes  when  the  local  system  is  considered. 

Under  more  realistic  conditions  (Condition  C)  the  flow  is  gener¬ 
ated  by  physical  forcing  as  well  as  truncation  error.  The  amount  of 
flow  generated  under  this  condition  with  the  standard  terrain  following 
coordinate  system  is  not  appreciably  different  (with  the  exception  of 
grid  point  5)  from  the  amount  generated  solely  by  truncation  error 
(Conditions  A  and  B).  Therefore,  under  realistic  conditions  it  may  be 
difficult  to  determine  variations  in  the  flow  field  due  solely  to 
physical  forcing.  This  is  not  the  case  with  the  local  terrain  follow¬ 
ing  system.  It  appears  that  the  truncation  error  has  been  reduced  to 
such  a  degree  (with  the  exception  of  grid  point  5)  that  variations  in 
the  flow  field  due  solely  to  physical  forcing  clearly  dominate  any 
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Figure  8.  Variation  of  the  maximum  U  component  with  time  above  grid 

point  5  for  the  local  {.)  and  standard  (x)  terrain  following 
coordinate  systems.  Left,  Condition  A  and  right.  Condition 


Variation  of  the  maximum  U  component  with  time  above  grid 
point  5  for  the  local  (.)  and  standard  (x)  terrain  follow 
coordinate  systems.  Condition  C. 


Figure  10.  Variation  of  the  maximum  U  component  with  time  above  grid 

point  14  for  the  local  (.)  and  standard  (x)  terrain  following 
coordinate  systems.  Left,  Condition  A  and  right.  Condition  B 


Variation  of  the  maximum  U  component  with  time  above  grid 
point  14  for  the  local  (.)  and  standard  (x)  terrain 
following  coordinate  systems.  Condition  C. 


Figure  12.  Variation  of  the  maximum  U  component  with  time  above  grid 
point  15  for  the  local  (.)  and  standard  (x)  terrain 
following  coordinate  systems.  Left,  Condition  A  and 
right.  Condition  B. 


Variation  of  the  maximum  U  component  with  time  above  grid 
point  15  for  the  local  (.)  and  standard  (x)  terrain 
following  coordinate  systems.  Condition  C. 


Figure  14.  Variation  of  the  maximum  U  component  with  time  above  grid 
point  18  for  the  local  (.)  and  standard  (x)  terrain 
following  coordinate  systems.  Left,  Condition  A  and 
right.  Condition  B. 


of  the  maximum  U  component  with  time  above  grid 
or  the  local  (.)  and  standard  (x)  terrain 
coordinate  systems.  Condition  C. 


Figure  16.  Variation  of  the  maximum  U  component  with  time  above  grid 
point  25  for  the  local  (.)  and  standard  (x)  terrain 
following  coordinate  systems.  Left,  Condition  A  and 
right.  Condition  B. 


Variation  of  the  maximum  U  component  with  time  above  grid 
point  25  for  the  local  (.)  and  standard  (x)  terrain 
following  coordinate  systems.  Condition  C. 


Variation  of  the  maximum  U  component  with  time  above  grid 
point  27  for  the  local  (.)  and  standard  (x)  terrain 
following  coordinate  systems.  Condition  C. 


Figure  22.  Variation  of  the  vertical  velocity  at  the  top  of  the  model 

with  time  above  grid  point  5  for  the  local  terrain  following 
coordinate  system  (.),  left,  and  for  the  standard  terrain 
following  coordinate  system  (x),  right.  Condition  C. 
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Figure  24.  Variation  of  the  vertical  velocity  at  the 
with  time  above  grid  point  14  for  the  loca 
coordinate  system  (.),  left,  and  for  the  s 
following  coordinate  system  (x),  right.  C 


standard  ter 
Condition  A. 
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me  above  grid  point  15  for  the  local  terrain  folloi 
late  system  (.),  left,  and  for  the  standard  terrain 
ng  coordinate  system  (x),  right.  Condition  C. 


m  of  the  vertical  velocity  at  the  top  of  the  model 
le  above  grid  point  18  for  the  local  (.)  and  standard 
ain  following  coordinate  systems.  Condition  C. 


Variation  of  the  vertical  velocity  at  the  top  of  the  model 
with  time  above  grid  point  27  for  the  local  (.)  and  standard 
(x)  terrain  following  coordinate  systems.  Condition  B. 


CONDITION  C  CONDITION  B  CONDITION 


Table  4.  A  coordinate  system  comparison  of  averages  of  |U  j 
from  selected  grid  points.  max 


Table  5.  A  coordinate  system  comparison  of  averages  of  |W.  j 
from  selected  grid  points. 


flow  associated  with  truncation  error. 

A  numerical  instability  was  apparently  caused  when  the  standard 
terrain  following  coordinate  system  was  employed.  The  presence  of  this 
instability  may  indicate  that  the  present  terrain  configuration  cannot 
be  resolved  by  the  standard  system  without  additional  smoothing.  If 
this  is  indeed  the  case,  the  local  system  enjoys  an  additional  advantage 
over  the  standard  system.  That  advantage  is  the  ability  to  resolve  a 
coarser  terrain  configuration  and  still  present  stable  solutions. 

Figures  36  through  47  present  an  opportunity  to  examine  the  hori¬ 
zontal  variation  of  the  truncation  error  produced  by  the  two  coordinate 
systems.  These  figures  show  the  integrated  vertical  velocity  field  at 
the  top  of  the  model  for  each  case  tested  after  integrations  of  eight 
and  eighteen  hours.  If  the  air  was  rising  the  barb  points  up;  if  the 

air  was  sinking,  the  barb  points  down.  The  plotting  convention  is  as 

85 

follows:  j"  =  lmm/sec,  f  =  55mm/sec,  f=  lOmm/sec,  and  j  =  85mm/sec 
(too  large  to  plot). 

7.2  Latitude 

7.2.1  Purpose 

The  purpose  of  this  test  is  to  determine  the  sensitivity  of  the 
forecast  vertical  motion  field  to  latitude. 

7.2.2  Test  Description 

As  discussed  in  Chapter  1,  the  phenomena  believed  responsible  for 
boundary  layer  convergence  are  the  diurnal  oscillation  of  eddy  stresses 
and  buoyancy  forces  above  sloping  terrain.  Closely  associated  with 
these  diurnal  forcings  is  the  inertial  oscillation  of  the  ageostrophic 
component  of  the  wind.  This  oscillation  is  linked  to  the  formation  of 


Figure  36.  Vertical  velocity  at  the  top  of  the  model  after  eight  hours 
with  the  local  terrain  following  coordinate  system.  Condition 
A.  The  plotting  convention  is  as  follows:  T  represents  rising 
motion,  1  represent'  sinking  motion,  f  =  lmm/sec,  r=  5mm/sec, 
f =  lOmm/sec,  and  j  =  85mm/sec  (too  large  to  plot). 


Figure  37.  Vertical  velocity  at  the  top  of  the  model  after  eight  hours 
with  the  standard  terrain  following  coordinate  system. 
Condition  A.  See  Figure  36  for  plotting  convention. 


Figure  39.  Vertical  velocity  at  the  top  of  the  model  after  18  hours 
with  the  standard  terrain  following  coordinate  system, 
Condition  A.  See  Figure  36  for  plotting  convention. 
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Figure  42.  Vertical  velocity  at  the  top  of  the  model  after  18  hours 
with  the  local  terrain  following  coordinate  system. 
Condition  B.  See  Figure  36  for  plotting  convention. 
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Figure  43.  Vertical  velocity  at  the  top  of  the  model  after  18  hours 
with  the  standard  terrain  following  coordinate  systems. 
Condition  B.  See  Figure  36  for  plotting  convention. 


4.5 


Vertical  velocity  at  the  top  of  the  model  after  eight  hours 
with  the  local  terrain  following  coordinate  system. 
Condition  C.  See  Figure  36  for  plotting  convention. 


Vertical  velocity  at  the  top  of  the  model  after  eight  hours 
with  the  standard  terrain  following  coordinate  system. 
Condition  C.  See  Figure  36  for  plotting  convention. 


Vertical  velocity  at  the  top  of  the  model  after  18  hours 
with  the  local  terrain  following  coordinate  system, 
Condition  C.  See  Figure  36  for  plotting  convention. 


Vertical  velocity  at  the  top  of  the  model  after  18  hours 
with  the  standard  terrain  following  coordinate  system. 
Condition  C.  See  Figure  36  for  plotting  convention. 
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low  level  jets  and  boundary  layer  convergence.  The  period  of  this 
oscillation  is  a  function  of  latitude  through  the  coriolis  parameter. 

Period  *  2ir/f  where  f  =  Coriolis  parameter 

f  =  2n  sin$  where  n  =  earth's  angular  velocity 

$  =  latitude 

The  lower  the  latitude,  the  longer  the  period;  therefore,  boundary  layer 
convergence  fields  formed  by  the  above  type  of  forcing,  would  tend  to 
occur  later  in  the  day  at  lower  latitudes  due  to  the  longer  period  of 
the  natural  oscillation.  It  should  also  be  mentioned  that  the  diurnal 
forcing  may  be  in  resonance  with  the  natural  period  of  oscillation  at 
certain  latitudes,  thereby  affecting  the  amplitude  of  the  winds  and 
possibly  the  boundary  layer  convergence. 

It  would  therefore  be  of  interest  to  examine  the  amplitude  and 
phasing  of  the  vertical  motion  fields  forecast  by  the  boundary  layer 
model  for  different  latitudes. 

Latitudes  38°  (f  =  .0000898),  30°  (f  =  .00007292),  and  22°  (f  = 
.00005463)  were  selected  for  the  tests.  The  initial  atmosphere  for  the 
model  contained  moisture  and  a  southerly  geostrophic  wind  of  lOm/sec. 

The  U  and  W  components  of  momentum  were  zero.  The  V  component  increased 
logarithmically  from  zero  at  the  roughness  height  to  lOm/sec  at  50m  and 
was  then  constant  with  height  to  the  top  of  the  model.  The  complete 
treatments  for  radiation  and  friction  were  used  in  these  tests.  The 
boundary  layer  model  was  executed  for  a  60  hour  forecast  period  at  lati¬ 
tudes  38°,  30°,  and  22°. 

7.2.3  Results 

The  amplitude  and  phase  of  the  maximum  vertical  motion  during  the 
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forecast  period  above  each  grid  point  is  depicted  in  Figures  48  through 
50.  The  plotting  convention  is  as  follows: 

A  vector  directed  from  the  north  denotes  a  midnight  maximum, 
while  one  directed  from  the  east  denotes  a  6  A.M.  (local 
time)  maximum,  etc.  Each  barb  (r)  denotes  lmm/sec;  each 
full  barb  (f)  denotes  5mm/sec.  Example:  The  amplitude 

and  phasing  for  the  example  is  1800  local  time  and  6mm/sec, 
respectively. 

Changes  in  phase  and  amplitude  are  noted  at  many  grid  points. 

A  closer  examination  is  made  of  five  grid  points;  the  results,  of 
which,  are  shown  in  Figures  51  through  55.  These  figures  contain  two 
parts.  The  upper  part  is  a  dial  which  depicts  the  phase  and  amplitude 
of  Wmax  for  each  latitude.  The  position  of  the  latitude  around  the 
dial  denotes  the  time  of  the  maximum  W  for  that  latitude  while  the 
distance  from  the  center  of  the  dial  denotes  the  amplitude.  The  ampli¬ 
tudes  are  plotted  in  mm/sec  with  the  scale  noted  on  the  dial.  The 
scale  changes  from  dial  to  dial  to  better  depict  the  particular  combin¬ 
ation  of  values  to  be  plotted.  The  lower  part  of  the  figure  shows  the 
relative  position  of  the  grid  point  in  question. 

Figures  51  through  55  show  a  change  of  phase  of  approximately  six 
hours  between  latitudes  22°  and  38°.  The  time  of  the  maximum  increases 
toward  lower  latitudes,  which  is  in  agreement  with  the  above  argument 
and  Paegle's  (1978)  linearized  solutions.  There  appears  to  be  no 
significant  change  in  amplitude  between  latitudes. 
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Figure  48.  Amplitude  and  phase  of  the  maximum  vertical  velocity  forecast  for 
38°  latitude.  The  plotting  convention  is  as  follows:  A  vector 
directed  from  the  north  denotes  a  midnight  maximum,  while  one 
directed  from  the  east  denotes  a  6  A.M.  (local  time)  maximum,  etc 
f=  lmm/sec  and  f*  5mm/sec. 


Amplitude  and  phase  of  the  maximum  vertical  velocity  forecast 
for  30°  latitude.  See  Figure  48  for  plotting  convention. 


Amplitude  and  phase  of  the  maximum  vertical  velocity  forecast 
for  22°  latitude.  See  Figure  48  for  plotting  convention. 
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Figure  51.  Amplitude  and  phase  of  the  maximum  vertical  velocity  fore¬ 
cast  at  grid  point  16  for  22°,  30°,  and  38°  latitude.  The 
position  of  the  latitude  around  the  dial  denotes  the  time 
of  the  max  W  for  that  latitude  while  the  distance  from  the 
center  of  the  dial  denotes  the  amplitude  in  mm/sec. 
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Figure  53.  Amplitude  and  phase  of  the  maximum  vertical  velocity  fore 
cast  at  grid  point  20  for  22°,  30°,  and  38°  latitude.  Se» 
Figure  51  for  plotting  convention. 


TERRAIN 

Amplitude  and  phase  of  the  maximum  vertical  velocity  fore¬ 
cast  at  grid  point  22  for  22°,  30°,  and  38°  latitude.  See 
Figure  51  for  plotting  convention. 


TERRAINS 


Figure  55.  Amplitude  and  phase  of  the  maximum  vertical  velocity  fore' 
cast  at  grid  point  25  for  22° »  30°,  and  38°  latitude.  Se< 
Figure  51  for  plotting  convention. 
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7.3  Radiation 

7.3.1  Purpose 

The  purpose  of  the  radiation  tests  is  to  determine  an  efficient 
and  economical  method  of  computing  the  radiation  fields. 

7.3.2  Test  Description 

Approximately  30  minutes  of  Univac  1108  CPU  time  is  required  for 
the  complete  boundary  layer  model  to  produce  a  four  hour  forecast  for 
only  seven  grid  points.  Radiation  calculations  were  completed  each 
time  step  and  accounted  for  much  of  the  CPU  time.  In  order  to  reduce 
the  computer  time  to  an  acceptable  level  without  adversely  affecting  the 
accuracy  of  the  model,  numerous  test  were  conducted.  These  tests  were 
to  check  the  effect  of  reducing  the  depth  of  the  radiational  model  as 
well  as  the  frequency  with  which  the  radiation  calculations  are  performed. 

7.3.3  Results 

The  top  of  the  radiational  model  was  lowered  from  nine  to  six 
kilometers.  This  change  resulted  in  the  elimination  of  over  30  compu¬ 
tational  levels.  The  vertical  grid  spacing  was  changed  from  100m  to 
86m  in  order  to  make  the  radiational  and  boundary  layer  models  more  com¬ 
patible.  Table  6  shows  the  effect  these  modifications  had  on  the  values 
of  some  of  the  more  important  radiational  parameters. 

The  next  step  was  to  test  the  sensitivity  to  the  frequency  with 
which  the  radiation  calculations  were  performed.  Numerous  tests  were 
conducted  in  which  the  radiation  calculations  were  done  as  frequently 
as  every  time  step  and  as  infrequently  as  every  four  hours.  A  test  was 
also  conducted  in  which  the  flux  divergence  was  held  constant. 

The  tests  in  which  the  radiation  is  calculated  every  two,  three, 
and  four  hours  used  the  complete  boundary  layer  model  as  described  in 
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Table  6.  Comparison  of  parameter  values  computed  from  radiational 
models  of  different  depths  and  vertical  grid  spacing 
(McClatchey  mid-latitude  summer  data  was  used  for  the 
calculations). 


Vertical  Grid 
Separation  (m) 


Model  Depth  (km) 


.01  m 
49.00  m 
306.00  m 
653.00  m 
1000.00  m 


100 

86 

9 

5 

3 

5 

2.424 

2.395 

2.309 

2.395 

31.61 

31.68 

31.98 

31.69 

196.1 

196.2 

196.4 

196.2 

164.49 

164.52 

164.42 

164.51 

-3.53 

-3.57 

-3.59 

-3.55 

-3.07 

-3.08 

-3.13 

-3.08 

-2.45 

-2.46 

-2.52 

-2.46 

-2.36 

-2.37 

-2.46  ; 

-2.37 

-2.34 

-2.37 

-2.46 

-2.37 
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Section  7.2  over  a  60  hour  forecast  period.  The  constant  flux  test 
was  also  conducted  using  this  model  with  the  following  specifications: 
(1)  global  radiation  updated  every  two  hours,  (2)  Q=0,  (3)  1^  *  3.4, 
and  (4)  F=0.  This  test  was  also  conducted  over  a  60  hour  forecast 
period.  The  effect  of  these  tests  upon  the  amplitude  and  phase  of  the 
vertical  motion  and  subsidence  fields  is  shown  in  Figures  56  through 
r3.  The  plotting  convention  is  described  in  Section  7.2.  Differences 
in  phase  and  amplitude  are  found  between  all  tests  with  the  greatest 
similarity  between  the  subsidence  fields  of  the  tests  where  radiative 
fluxes  were  calculated  every  three  hours  and  every  four  hours. 

The  thermal  structure  was  closely  examined  through  an  analysis  of 
the  diurnal  temperature  variations  at  grid  points  15,  18,  and  22  (see 
Figure  64).  Temperatures  at  the  roughness  height,  49m,  133m,  and 
480m  were  plotted  over  a  48  hour  period  as  forecast  in  the  two,  three, 
and  four  hour  radiational  tests  (Figures  65  through  76).  The  plotting 
convention  is  as  follows: 

+  Radiation  update  at  two  hour  intervals 
0  Radiation  update  at  three  hour  intervals 
a  Radiation  update  at  four  hour  intervals 

Temperature  curves  from  the  above  radiational  tests  are  similar  at 
each  level . 

The  tests  described  in  this  section  did  not  point  to  a  "best"  fre¬ 
quency  for  the  radiational  calculations;  therefore  a  two  hour  frequency 
was  selected  since  this  was  the  highest  frequency  considered  econom¬ 
ically  feasible.  This  frequency  adequately  resolves  temperature 
oscillations  with  periods  greater  than  four  hours.  Based  on  extensive 
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Amplitude  and  phase  of  the  maximum  vertical  velocity  witl 
radiation  fields  updated  every  two  hours.  See  Figure  48 
for  plotting  convention. 
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Amplitude  and  phase  of  the  maximum  vertical  velocity  wi 
radiation  fields  updated  every  three  hours.  See  Figure 
48  for  plotting  convention. 


Amplitude  and  phase  of  the  maximum  vertical  velocity  with 
radiation  fields  updated  every  four  hours.  See  Figure  48 
for  plotting  convention. 


Figure  59.  Amplitude  and  phase  of  the  maximum  vertical  velocity  with 
infrared  radiation  a  constant.  See  Figure  48  for  plotting 
convention. 


Figure  60.  Amplitude  and  phase  of  the  maximum  subsidence  with 

radiational  fields  updated  every  two  hours.  See  Figure  48 
for  plotting  convention. 


Figure  61.  Amplitude  and  phase  of  the  maximum  subsidence  with 
radiational  fields  updated  every  three  hours.  See 
Figure  48  for  plotting  convention. 


Figure  62.  Amplitude.  and  phase  of  the  maximum  subsidence  with 
radiation  fields  updated  every  four  hours.  See 
Figure  48  ior  plotting  convention. 


Figure  63.  Amplitude  and  phase  of  the  maximum  subsidence  with 
infrared  radiation  a  constant.  See  Figure  48  for 
plotting  convention. 


Figure  64.  Grid  points  selected  for  a  study  of  the  radiational 
effects  on  the  thermal  structure. 
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Figure  65.  Diurnal  temperature  curves  from  grid  point  15.  Height  is 
.01m  .  +  =  radiation  update  at  two  hour  intervals,  •  = 
radiation  update  at  three  hour  intervals  and  A  =  radiation 
update  at  four  hour  intervals. 
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LOCAL  TIME  (HRS) 

Figure  66.  Diurnal  temperature  curves  from  grid  point  15.  Height  is 
49m.  See  Figure  65  for  plotting  convention. 


Diurnal  temperature  curves  from  grid  point  15.  Height  is 
133m.  See  Figure  65  for  plotting  convention. 


Diurnal  temperature  curves  from  grid  point  15.  Height  is 
480m.  See  Figure  65  for  plotting  convention. 


Diurnal  temperature  curves  from  grid  point  IS.  Height  i 
•01m.  See  Figure  65  for  plotting  convention. 


Diurnal  temperature  curves  from  grid  point  18.  Height  is 
49m.  See  Figure  65  for  plotting  convention. 


Diurnal  temperature  curves  from  grid  point  18.  Height  is 
133m.  See  Figure  65  for  plotting  convention. 


Diurnal  temperature  curves  from  grid  point  18.  Height  is 
480m.  See  Figure  65  for  plotting  convention. 


Diurnal  temperature  curves  from  grid  point  22.  Height  is 
•01m.  See  Figure  65  for  plotting  convention. 


Diurnal  temperature  curve  from  grid  point  22.  Height  is 
133m.  See  Figure  65  for  plotting  convention. 


Diurnal  temperature  curves  from  grid  point  22.  Height  is 
480m.  See  Figure  65  for  plotting  convention. 


radiational  tests  with  the  model,  it  is  felt  that  this  limitation  is 
not  an  area  of  concern. 

7.4  Soil  Type 

7.4.1  Purpose 

The  purpose  of  this  test  is  to  determine  the  sensitivity  of  the 
boundary  layer  vertical  motion  field  to  soil  type. 

7.4.2  Test  Description 

Linearized  results  of  Paegle  (1978)  indicated  a  sensitivity  of  the 
boundary  layer  convergence  fields  to  the  evolution  of  the  low-level 
temperature  field.  The  question  arises  as  to  the  effect  of  soil  type 
on  the  interacting  low-level  temperature  field  and  momentum  boundary 
1 ayer . 

Three  very  different  soil  types  were  tested  using  the  complete 
boundary  layer  model  as  described  in  Section  7.2.  The  amplitude  and 
phase  of  the  vertical  motion  field  as  well  as  the  diurnal  temperature 
patterns  from  several  grid  points  were  analysed. 

7.4.3  Results 

The  soils  were  classified  by  conductive  capacity  which  is 
This  term  is  a  measure  of  the  soil's  efficiency  in  diffusing  heat 
away  from  the  surface  and  retaining  it,  (the  larger  the  capacity,  the 
more  efficient  the  soil).  The  three  soils  tested  and  their  identifying 
parameters  are  listed  in  Table  7. 

The  amplitude  and  phase  of  the  vertical  velocity  fields  as  a  func¬ 
tion  of  the  soil  parameters  are  shown  in  Figures  77,  78,  and  79.  A 
summary  of  this  data  is  contained  in  Table  8.  Column  A  of  Table  8 
lists  the  soil  types  that  are  being  changed;  Column  B,  the  average 
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Table  7.  Parameters  for  the  three  soils  tested. 


Soil  Type 

- — 

pscs(cal  m”3oK) 

Ks(m^sec"^ ) 

Conductive  Capacity 

Dry  Sand 

3.0  x  105 

1.3  x  10-7 

108 

Clayland 

5.6  x  105 

1.2  x  10'7 

194 

Sandy  Clay 

5.9  x  105 

3.7  x  10‘7 

357 

Table  8.  Effects  of  altering  soil  type  on  amplitude  and  phase  of 
boundary  layer  vertical  velocity  fields 


Change  in 

Soil  Type 

Average  Difference  In 
Phase(hrs)/Number 
of  Grid  Points 

Average  Difference  In 
Ampl itude(mm/sec)/ 
Number  of  Grid  Points 

CP  to  DS 

3/28 

-9/3 

8/30 

-3/1 

SC  to  CP 

4/15 

-7/7 

6/31 

— 

SC  to  DS 

6/21 

_  _ i 

14/31 

DS  =  Dry  Sand,  CP  =  Clayland  Pasture;  SC  =  Sandy  Clay 

increase  and  decrease  in  the  phase  of  the  maximum  vertical  velocity  and 
the  grid  points  involved  as  a  result  of  the  soil  change;  and  Column  C, 
the  same  as  Column  B  except  this  data  pertains  to  the  amplitude  of  the 
maximum  vertical  velocity. 

A  comparison  of  the  amplitudes  and  phases  of  the  maximum  vertical 
velocity  fields  for  different  soil  types  is  made  at  selected  grid 
points  (Figures  80  through  87).  The  plotting  convention  is  explained 


Figure  77.  Amplitude  and  phase  of  maximum  vertical  velocity  with  dry 
sand  soil  (conductive  capacity  =  108).  See  Figure  48  for 
plotting  convention. 


Figure  78.  Amplitude  and  phase  of  maximum  vertical  velocity  with 
clayland  pasture  (conductive  capacity  =  194).  See 
Figure  48  for  plotting  convention. 


Figure  79.  Amplitude  and  phase  of  maximum  vertical  velocity  with 
sandy  clay  (conductive  capacity  =  357).  See  Figure 
48  for  plotting  convention. 
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Figure  80.  Amplitude  and  phase  of  maximum  vertical  velocity  for  sandy 
clay  (X),  clayland  pasture  (P),  and  dry  sand  (S)  for  grid 
point  4.  See  Figure  51  for  plotting  convention. 
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Figure  81.  Amplitude  and  phase  of  maximum  vertical  velocity  for  sandy 
clay  (X),  clayland  pasture  (P),  and  dry  sand  (S)  for  grid 
point  8.  See  Figure  51  for  plotting  convention. 
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Amplitude  and  phase  of  maximum  vertical  velocity  for  sandy 
clay  (X),  clayland  pasture  (P),  and  dry  sand  (S)  for  grid 
point  15.  See  Figure  51  for  plotting  convention. 


■cm  »i  n:annirtli 


•  i-irtruaaiwi 


i i  m«wmm  nmumsHM 


TERRAIN 


Figure  84.  Amplitude  and  phase  of  maximum  vertical  velocity  for  sand 
clay  (X),  clayland  pasture  (P),  and  dry  sand  (S)  for  grid 
point  20.  See  Figure  51  for  plotting  convention. 


Figure  85.  Amplitude  and  phase  of  maximum  vertical  velocity  for  sandy 
clay  (X),  clayland  pasture  (P),  and  dry  sand  (S)  for  grid 
point  22.  See  Figure  51  for  plotting  convention. 


TERRAIN 

87.  ...plitude  and  phase  of  maximum  vertical  velocity  for  sandy 
clay  (X),  clayland  pasture  (P),  and  dry  sand  (S)  for  grid 
point  28.  See  Figure  51  for  plotting  convention. 
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in  Section  7.2.  These  figures  and  Table  8  clearly  indicate  that  a 
change  in  soil  type,  such  that  the  conductive  capacity  is  altered,  may 
result  in  a  substantial  change  in  the  amplitude  and  phase  of  the  maxi¬ 
mum  vertical  velocity.  In  the  extreme  case  of  changing  from  a  sandy 
clay  to  a  dry  sand  soil,  the  maximum  amplitude  occurred  an  average  of 
six  hours  later  for  21  of  the  31  grid  points  and  the  amplitude  in¬ 
creased  an  average  of  14mm/sec  at  each  grid  point. 

Figures  88  through  99  show  the  relationship  between  soil  type  and 
diurnal  temperature  pattern.  The  following  plot  convention  was  used: 
(1)  •  clayland  pasture,  (2)  +  dry  sand,  and  (3)  A  sandy  clay.  The 
larger  the  conductive  capacity  of  the  soil  the  less  variation  in  the 
temperature  field.  These  figures  also  indicate  that  a  change  in  soil 
type  does  not  only  affect  the  very  low  levels  of  the  boundary  layer  but 
also  at  least  above  500m  where  a  four  to  five  degree  temperature  dif¬ 
ference  is  noted  between  soil  types. 
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Figure  88.  Diurnal  temperature  curves  for  grid  point  15.  Height  is 
See  Figure  65  for  plotting  convention. 


Diurnal  temperature  curves  for  grid  point  15 
See  Figure  65  for  plotting  convention. 


Figure  90.  Diurnal  temperature  curves  for  grid  point  15.  Height  is  133m. 
See  Figure  65  for  plotting  convention. 
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Diurnal  temperature  curves  for  grid  point  15.  Height  is  480m. 
See  Figure  65  for  plotting  convention. 


Diurnal  temperature  curves  for  grid  point  18.  Height  is 
See  Figure  65  for  plotting  cuivention. 


Diurnal  temperature  curves  for  grid  point  18. 
See  Figure  65  for  plotting  convention. 


Diurnal  temperature  curves  for  grid  point  18.  Height  is  133m. 
See  Figure  65  for  plotting  convention. 


Diurnal  temperature  curves  for  grid  point  18.  Height  is  480m. 
See  Figure  65  for  plotting  convention. 


Diurnal  temperature  curves  for  grid  point  22.  Height  is 
See  Figure  65  for  plotting  convention. 
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Figure  97.  Diurnal  temperature  curves  for  grid  point  22.  Height  is  49m. 
See  Figure  65  for  plotting  convention. 
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Figure  98.  Diurnal  temperature  curves  for  grid  point  22.  Height  is  133m. 
See  Figure  65  for  plotting  convention. 
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Diurnal  temperature  curves  for  grid  point  22.  Height  is  480m 
See  Moure  65  for  plotting  convention. 


CHAPTER  8 


CONCLUSIONS 

The  intent  of  this  research  has  been  to  develop  a  boundary  layer 
model  with  sufficient  resolution  and  physics  to  define  the  interacting 
thermal  and  momentum  boundary  layers  on  diurnal  time  scales  above 
synoptic  scale  sloping  terrain  and  thermally  inhomogeneous  surfaces. 

It  is  felt  that  this  model  will  add  understanding  to  the  physical 
processes  and  simulation  requirements  of  diurnal  boundary  layer  motions. 

Proper  resolution  of  the  terrain  configuration  is  essential  in 
models  of  this  type  since  the  terrain  slope  is  one  of  the  primary 
forcing  mechanisms  for  the  diurnal  boundary  layer  motions.  This  re¬ 
search  made  significant  progress  in  that  area.  A  local  terrain 
following  coordinate  system  was  developed  which  has  characteristics 
superior  to  the  terrain  following  coordinate  systems  commonly  used. 

The  local  coordinate  system  produces  truncation  errors  which  are 
significantly  less  than  those  produced  by  the  standard  systems. 

Variations  in  the  flow  field  due  to  physical  forcing  are  discern- 
able  from  those  variations  due  to  truncation  error  in  the  local 
coordinate  system.  This  is  not  always  the  case  with  the  standard 
coordinate  systems.  In  the  latter,  it  is  not  always  possible  to  deter¬ 
mine  that  variations  in  the  flow  field  are  caused  by  physical  forcing 
or  truncation  error.  The  local  coordinate  system  enjoys  an  additional 
advantage  over  the  standard  systems  in  that  it  can  apparently  resolve 
a  coarser  terrain  configuration  and  still  provide  stable  solutions. 
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The  sensitivity  of  the  complete  model  was  tested  against  para¬ 
meters  which  had  shown  a  certain  amount  of  sensitivity  in  earlier,  less 
complete  models.  The  phenomena  believed  responsible  for  boundary  layer 
convergence  are  the  diurnal  oscillations  of  eddy  stresses  and  buoyancy 
forces  above  sloping  terrain.  The  argument  can  be  made  that  if  bound¬ 
ary  layer  convergence  fields  are  formed  by  the  above  type  of  forcing, 
their  formation  should  occur  later  in  the  day  at  lower  latitudes. 

Also  the  diurnal  forcing  may  be  in  resonance  with  the  natural  period 
of  oscillation  at  certain  latitudes,  thereby  affecting  the  amplitude 
of  the  winds  and  possibly  the  boundary  layer  convergence.  Forecasts 
by  the  complete  model  did  reveal  a  sensitivity  to  latitude.  The 
phase  of  the  maximum  vertical  motion  occurred  six  hours  later  between 
38°N  and  22°N.  These  results  are  in  agreement  with  Paegle's  (1978) 
linearized  solutions.  The  model  did  not  indicate  a  significant  change 
in  amplitude  between  latitudes. 

Radiation  calculations  in  the  complete  model  were  very  time  con¬ 
suming.  Many  sensitivity  tests  were  conducted  to  find  ways  to  reduce 
this  time.  These  tests  revealed  that  a  radiational  model  5  KM  thick 
was  as  effective  in  calculating  boundary  layer  global  radiation  and 
flux  divergence  as  a  model  9  KM  thick. 

The  boundary  layer  model  is  apparently  not  very  sensitive  to  flux 
divergence  calculations;  in  fact,  tests  conducted  to  determine  the  best 
frequency  for  radiation  calculations  did  not  show  significant  differ¬ 
ences  between  every  two,  three,  and  four  hours.  A  two  hour  frequency 
was  selected  since  it  was  the  highest  frequency  considered  ecomonically 
feasible.  This  frequency  adequately  resolves  temperature  oscillations 
with  periods  greater  than  four  hours  and  is  therefore  not  a  serious 
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limitation  for  the  model.  Caution  should  be  expressed  when  the  sensi¬ 
tivity  of  the  boundary  layer  model  to  flux  divergence  is  considered. 
Approximations  in  the  radiational  model  and  emissivity  values  resulted 
in  cooling  rates  which  were  less  than  those  calculated  in  more  complete 
radiational  models. 

Little  has  been  reported  in  the  literature  on  the  relationship  be¬ 
tween  soil  type  and  boundary  layer  convergence.  Sensitivity  tests  of 
the  complete  model  with  three  very  different  soil  types  shows  a  defi¬ 
nite  connection  between  the  amplitude  and  phase  of  the  maximum  vertical 
velocity  and  soil  type.  In  the  extreme  case  of  changing  from  a  sandy 
clay  to  a  dry  sand  soil,  the  phase  occurs  six  hours  later  and  the 
amplitude  increases  14  nm/sec.  It  is  apparent  from  these  results  that 
future  boundary  layer  models  should  seriously  consider  the  inclusion 
of  soil  type. 

The  present  results  tends  to  support  the  viewpoint  that  the  pro¬ 
nounced  diurnal  control  upon  air  mass  thunderstorms  over  the  United 
States  resides  in  boundary  layer  processes.  This  control  is  not  simply 
in  reponse  to  diurnal  low  level  stability  changes,  but  rather  through 
the  boundary  layer  convergence  field  which  is  significantly  determined 
by  terrain  slopes  and  curvatures  as  gentle  as  1/500  and  2/106  km-^ , 
respectively.  Such  terrains  are  found  over  the  central  plains  of  the 
United  States.  Here,  the  response  of  W  on  the  order  of  1  mm/sec  to 
10  mm/sec  at  1  KM  corresponds  to  layer  average  divergences  of  order 
10"6  sec-1  to  10'5  sec"^ .  These  magnitudes  are  similar  to  those  asso¬ 
ciated  with  synoptic  scale  weather  systems  and  suggest  that  precise 
timing  of  thunderstorm  activity  may  be  predicted  only  with  a  coupled 
boundary  layer-free  tropospheric  model. 
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In  the  past,  major  problems  in  the  execution  of  such  a  model  have 
involved  excessive  numerical  truncation  error  in  the  vicinity  of  slop¬ 
ing  terrain  and  uncertainty  with  respect  to  the  most  relevant  boundary 
layer  processes.  Both  of  these  problems  have  been  addressed  in  this 
research  with  some  success.  It  appears  that  the  next  major  obstacle 
is  the  availability  of  a  sufficiently  powerful  computer  to  process  a 
fully  three  dimensional  coupled  model.  Based  on  the  computer  require¬ 
ments  of  the  present  two  dimensional  model,  this  may  be  technologically 
feasible  on  the  current  largest  computers.  However,  it  appears  that 
the  most  profitable  approach  to  explore  the  scientific  questions  and 
to  develop  the  numerical  techniques  is  with  a  two  dimensional  version 
as  done  here. 


APPENDIX  A 


DERIVATION  OF  TERRAIN  FOLLOWING 
COORDINATE  SYSTEM 


In  order  to  derive  the  equations  it  is  necessary  to  connect  the 
coordinates  of  one  point  in  two  Cartesian  reference  frames  possessing 
the  same  origin.  The  equations  are  developed  for  the  general  case  in 
which  both  systems  are  oblique  and  then  for  the  special  case  in  which 
one  system  is  oblique  and  the  other,  orthogonal. 

The  systems  used  are  0-XYZ-aBy  and  O-X-iY^Z^-ct-jS-jT-j  with  the  units 
the  same  in  each  system  as  shown  in  Figure  100.  The  cosines  of  the 
angles  formed  by  the  axes  of  these  systems  are  indicated  in  Table  9. 


Table  9.  Table  of  Cosines 


0 

X 

Y 

Z 

X1 

£1 

ml 

nl 

Y1 

l2 

"2 

*3 

m3 

n3 

Example:  cos  <  YOZ^  = 

Consider  a  closed  broken  line  which  leads  from  0  to  P  along  the 
edges  0PX,  P^y,  and  pxyP  the  sYstem  0_XYZ  and  which  returns  from 
P  to  0  along  the  edges  PP^  Y  »  Px  y  PX-|  ’  px^°  °'fr  SYStem  0_xi  Yizi 


as 


* 


150 


shown  in  Figure  101.  The  closed  broken  line  is  projected,  in  turn,  on 
the  axes  OX,  OY,  OZ  and  then  on  the  axes  OX^ ,  OY^  ,  0Z-| .  The  following 
theorems  are  useful : 

Theorem  1.  The  projection  of  a  segment,  AB,  of  a  directed  line, 
£,  upon  a  directed  line,  m,  is  equal  to  the  product  of  AB  by  the  co¬ 
sine  of  the  angle  between  the  two  directed  lines. 

Theorem  2.  The  sum  of  the  projections  upon  a  directed  line,  m, 
of  the  segments  of  a  closed  broken  line  in  space  is  equal  to  zero. 
Projecting  the  closed  broken  line  onto  OX,  OY,  and  OZ  yields  the  fol¬ 
lowing  system: 

X  +  Y  cosy  +  Z  COS6  -  X.|£-|  -  Y-j£2  -  Z-|£3  *  0 

X  cosy  +  Y  +  Z  cosa  -  X-jm-j  -  Y^m2  -  Z^m^  =  0 

X  cose  +  Y  cosa  +  Z  -  X^  -  Y^n2  -  Z^n^  =  0 

Projecting  the  closed  broken  line  onto  0X^ ,  OY^ ,  and  OZ^ ,  yields  the 

following  system: 


X£1 

+ 

Ym.j 

+ 

Z"1  -  X1 

-  Y 1  COSy^ 

-  zl 

COSB1 

=  0 

/  ?.2 

+ 

Ym2 

+ 

X 

1 

C\1 

C 

COSy-j  -  Y1 

-  zl 

COSa-j 

=  0 

XJI3 

+ 

Ym3 

+ 

Zn3  -  X] 

cose-|  -  Y-| 

COSa-j 

-  Z1 

=  0 

The  first  system  has  a  unique  solution  for  X  ,  Y,  Z  in  terms  of  X^ ,  Y1 , 
Z1  while  the  second  system  has  a  unique  solution  for  X-j ,  Yp  Z^  in 
terms  of  X,  Y,  Z.  Both  systems  are  solved  using  Cramer's  rule. 

The  derivation  up  to  this  point  follows  Dresden  (1930).  A  depar¬ 
ture  from  his  derivation  is  taken  when  the  following  specifications  are 
made: 

a.  Let  the  system  0-X^Y^Z^-a^e-|yi  be  orthogonal;  therefore  a^, 

8-j ,  y.|  are  right  angles. 


4 


Figure  101.  Closed  broken  line  which  goes  from  0  to  P  using  the  system 
0-XYZ  and  returns  from  P  to  0  by  way  of  the  system  O-XjY-.Z,. 
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b.  Let  the  system  O-XYZ-agy  be  our  "terrain  following  system" 
where  a  and  6  are  right  angles.  OX  rotates  in  the  X^Z-j  plane  and  is 
perpendicular  to  OY^.  OY  rotates  in  the  Y-j Z-j  plane  and  is  perpendicu¬ 
lar  to  OX-j  as  seen  in  Figure  102.  After  making  these  specifications, 
cosa-j  =  cos^i  =  cosy^  =  cosg  =  cosa  =  0.  Our  new  table  of  cosines 
appears  in  Table  10. 


Table  10.  Table  of  Cosines 


0 

X 

Y 

Z 

X1 

£1 

m-|=0 

nl 

Y1 

£2=0 

nu 

n2 

Z1 

n3 

The  two  systems  of  equations  are: 


Y1  =  m2Y  +  n2Z 

Z.|  =  «.3X  +  m^Y  +  n3Z 

£•,  m9  cosy  ju  -  m,  cosy 

X  =  ( - L-)X  -  (-2 - — )Y,  +  (-5 - - )Z, 

1-COS  y  1-COS  y  1-COS  y 

-Z-,  cosy  m9  m,  -  i~  cosy 

Y  =  (  p  )X,  +  ( - %-)Y,  +  (-* - - >Z1 

1-COS  y  1-COS  y  1-COS  y 

Z  ■  n1X1  +  n2Y1  +  n^ 


The  following  definitions  are  made  (refer  to  Figure  103  for  clarity): 
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£.|  =  COS  <  XOX-j  E  cos  ® 
n.|  =  cos  <  ZOX-j  =  cos  e 
m2  =  cos  <  Y0Y.J  s  cos  p 
n2  =  cos  <  ZOY^  i  cos  u 


z 2  =  cos  <  X0Z-|  =  sin  ® 

m^  =  cos  <  Y0Z-|  e  sin  p 

=  cos  <  ZOZ^  e  cos  x 


It  is  necessary  to  express  y,  x,  e,  and  y  in  terms  of  ®  and 
The  following  relations  can  be  established  (e.g.  Snyder,  1914): 


cosy  =  sin®  si n <}> 
1 


cosx  = 


cose  = 


^tan^®  +  sec^tp 
-tan® 

J  2  2 

tan  ffl  +  sec  4> 


COSy  = 


-tan<t> 


tan^®  +  sec% 


The  two  systems  of  equations  written  in  terms  of  only  ®  and  <t>  are  now: 


X-j  =  coseX  - 


tan® 


J  2  2 

tan  ®  +  sec 


Y1  =  cos<t>Y 


tan<f> 


j  5  2~ 

tan  ®  +  sec  <j> 


Z-j  =  sin®X  +  sin$Y  + 


J — 2 - 2 

tan  ®  +  sec  4 


X  -  [ 


2  2 

1-sin  ®sin  4 


lY  rcos4sin®sin4>n„  ,rsin®(1-sin  P)7  7 
JAl"L  p  o  JVL  p  2“*i 
1  -si n  ®si n  1-sin  ®sin  4> 

Y  „  [-C05tein»sir'.t.]x  c - cosj,  ]y  ^sintn-slr^j  , 

l-sin^®sin%  1  -sin^®sin%  l-sin^fisinS  1 


Z  -  [ 


-tan® 


]Xn-[- 


tanj) 


]Y,+[- 


1 


j  2  ^  Ly  2  T~  J  ^  ~ J  2"  T 

tan  ®+sec  4  tan  ®+sec  4  tan  SH-sec  p 


1  h 


Using  the  chain  rule  : 


a(  )  -  j (  )  r  cos9  i  3(  )  r  tan9 

3X,  3X  1  c<«2ifci«2  r  / - 9 - 5 — 

1  "S1n  951  $  /tan2®fsec24, 

3(  )  -  3(  )  r  COSd,  1  3  (  )r  tana 

3V  IV  I-  5 5^J  i7  L  _ 


3Y,  3Y  1  J  3Z 

I  1  “S i n  ®s in  ^ 


'  •  3,M  ^,n  *  tan  ®<-sec % 

il  )  _  a(  )  rsinS(1-sin26kf  3(  )rsind>(l-sin28)-,  a( 

3Z,  3X  1  e-i„2*.  2  J  3Y1!  2  J  C 

I  l-sin  ©sin  (j)  1-sin  ©sin  <p 


Hx1‘l)t^cose*iz^si' 


iL  1  =  iLl  rn^A  +  3(  )  cin* 

9Y  3Y^  C0S*  +  3Z^  S1"* 

li  )  -  3(  ) r  tanfl  1  _ ^P 

3Z  3XT1  J - 5 - 5-  J  lYT^- 

tan^GH-sec^ 


tan^®+sec^<)> 


_I[ _ tanfr  -|+  3_(  )r_  _1_ 

1  J  o  "  p  3Z-i  /  p  5  *" 
tan  »*-sec  <j>  tan  AM-sec  $ 


The  formulation  of  expressions  for  ,  V-j ,  and  W-j  where 
dX,  dY,  dZ, 

ui  =  ~&t  •  vi =  -ar  *  and  wi  =  "dr 

requires  much  manipulation  which  is  omitted  here  for  brevity.  The 
results  are: 


U,  =  cosuu - isn!_  H  -(sinax  +  ■ . ..  z)  a 

(tanVsec2*)3'2  dt 


x.  /  tan®tan4>secS  „x  d<j> 

Z)  d‘ 


V1  =  cos<i>V  - 


tan  ®+sec% 


2  2 

,,  u  ,  sec  <j>sec  ®  ,x  d* 

“  ‘(S’n  V  (^Lc!,|«2  Z)  * 


,  tanijitanfisec^®  -,x  d® 
+  - 7 - 5—372-  Z)  dt 


(tan  fH-sec  <j>) 
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After  much  manipulation  the  equations  are  formed  and  then  scaled 
according  to  the  following  factors: 

a.  Time  Scale  (T)  -  The  diurnal  time  scale  was  used;  therefore 
T~0(10^)  =  1/f  where  f  was  the  coriolis  parameter. 

b.  Depth  Scale  (D)  -  determined  by  the  physical  depth  of  the 

O 

model  which  was  1.5  km  '0(10  )  meters. 

c.  Horizontal  Wind  Speed  (V)  -  0(10)  meters  per  second. 
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d.  Length  Scale  (L)  -  determined  by  1/4  of  the  wavelength  of  the 
type  of  motion  which  could  be  resolved  on  the  model  grid.  For  a  grid, 
spacing  of  2°,  L~0(103)  meters. 

e.  Vertical  Scale  (W)  -  0(10  )  meters  per  second. 

f.  fx’  fy*  fx’  ff  "  <  This  was  based  on  a  smoothed  ter¬ 

rain. 

G.  e,<|>  -  0(5  10"3) 

The  pressure  gradient  term  requires  special  attention.  Using  the 
U  component  of  the  momentum  equation  as  an  example,  the  pressure  gradi¬ 
ent  term  is  linearized  about  a  basic  state  in  which  F  and  p  are  con¬ 
stant  in  X-j  and  t.  In  the  terrain  following  coordinate  system  we  have 
the  following: 


1  lE_  _  1 3ji  r 
P  ax1  "  -  1 


I  3P 

p  ax 


cos® 


2  2  J 
1-sin  ®sin  4> 


1  3P  r 
p  az  L 


tan® 


J ?  9 

tan  ©<-sec  < 


1  aF 

-  3X 


cos® 


1  -sin 


_ ]  _  1  r_ 

,in2,J  -  az  L y 

,m  4>  p  /. 


tan® 


2  2 
tan  ®+sec  < 


Note  that 


aP ' 

aZ 


■p  g  cosa  = 


JO_ 


/  2  2 

tan  ®+sec  < 


Therefore 


1  3P  ±  1  9P 1  r  COS®  -I  +  p!  „  r  tan®  -1 

p  aX^  -  oX  i-sin2©sin2i>  p  tan3®^sec24> 


The  equations  can  now  be  written  in  the  local  "terrain  following 
coordinates". 
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2  2 

3U  r  COS  ®  T  ,,  3U  ,  r  COS  *  n  „  all  ,  r  1 

Tf  +  L  p  p—J  U  rjr  +  L  r~J  V  ry  +L  p  “  p  j"ry 

l-sirr®sin*<|>  3A  1  -si  n^®s  i  n^4>  3T  tan^®+sec  $  3Z 


■]Wtt 


cos4> 

cos® 


fV  =  —  [- 


_ I _ .1  llPl.r 

;_2*_.._2,J  -  aX  L 


1  — s i n  ©sin  4>  p 


sin®  i  o' 

2 - 2 - 'J  -  9 

cos  Stan  0+1  p 


1 


3U, 


cos®  aZ 


lr+[ 


2 

cos  ® 


l-sin2Ssin2 


_i  u  —  +  r — -~-s-  - _ i  v  —  +r 

kJ  3X  L,  e.2.  .  Zj  V  *v 

Ji  1-sin  ®sin  0 


sY  tan2&+sec 


nw  31 
2jw  aZ 


+  COSi  fy  x  r. 
cos0  L 


1 


_ .1  1  ILL  _  r  ilHA  _  .1  £L_  a 

.._2^.._2;J  -  aY  Lr~2;_'  2,  j  -  3 


1-sin  Ssin  0  p 


cos  0tan  p 


+  — ! —  —  37) 

cos 0  az  n 


f?-+  [;— g|"--  p  ]  U  If  +  ]  V  |a  +  [-.--J..— -]w|9. 


3t  L1-sin2®sinV 


?  ? 

1  -si n  ©si n  4> 


aY 


?  2  J  ^7 

tan  ®+sec  <t  dL 


■  h  <Ke 


ae  +  r  cos  ® 
n  l-sin29sin2 


-]u|f+[ 


cos 


-3  v!t+  [ 


l-sin2®sin2 


30. 

p  p  J"^7 

tan  ®+sec  0 


-]W 


Q  +  A  (Ke  W 


[- 


1 


tan  ®+sec  $ 


_ pi  aw 

2  J  aZ 


=  -  [ 


2 

cos  fl 


2  2 
1-sin  ®sin  < 


i  au  r  cos  0  t  av 
J  3X  "  l-sin2®sin20  3Y 


The  two  dimensional  .  implified  version  of  these  equations  are  used  in 
the  model  as  explained  in  Chapter  3. 
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APPENDIX  B 


CALCULATING  SPECIFIC  HUMIDITY  FROM 
DEW  POINT  TEMPERATURE 

The  derivation  begins  with  the  Clausius-Clapeyron  equation 
(Rogers,  1976). 

des<T>  .  y 
dT  T(av-a£) 

where  e$  =  saturation  vapor  pressure 
L^v  =  latent  heat  of  vaporization 
av,a£  =  specific  volumn  for  vapor  and  liquid  phases  of  water 
Since  a  >>a  and  e  a  =R  T,  the  Clausius-Clapeyron  equation  can  be 

*  O  V  V 

written  in  approximate  form  as 

1  des  .  Ltv 
es  dT  RVT2 

where  Ry  =  individual  gas  constant  for  water  vapor. 

The  above  equation  is  integrated  from  some  reference  state  (ep, 
Tq)  to  a  final  state  (es,T).  Over  the  interval  of  integration,  L£y  is 
considered  constant  (L  =17  ).  Integration  yields 

es(T)  ;  e0  exp[^-  -  f)l 

Since  saturation  occurs  at  the  dew  point  (Tp),  the  above  equation  can 
be  interpreted  as 
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e  £  es^TD^  =  6,11  exp  ^ 273  "  T0>] 

where  the  reference  state  is  the  triple  point.  Appropriate  values  for 
r  and  Ry  are  from  Rogers  (1976)  and  Dutton  (1976) 

q  *  • 622(f) 


Substituting  the  above  expression  for  e  yields  the  final  result. 

„  .  r,,  ,6.11  nfnr  C597  -  -285(10  -  273)]  ,  , 

q  .622  {  p  exp[  ^273  " 


where  P  =  atmospheric  pressure  in  mb 
Td  =  dew  point  temperature  in  °K. 


REFERENCES 


Anthes,  R.  A.,  T.  T.  Warner  and  A.  L.  McNab,  1977:  Development  and 
testing  of  a  mesoscale  primitive  equation  model  at  Penn  State. 
Presented  at  the  Third  Conference  on  Numerical  Weather  Prediction 
of  the  A.  M.  S.,  Omaha,  Nebraska. 

Blackadar,  A.  K.,  1962:  Vertical  distribution  of  wind  and  turbulent 
exchange  in  a  neutral  atmosphere.  J.  Geophy.  Res.,  67,  3095- 
3102. 

Blackadar,  A.  K. ,  1957:  Boundary  layer  wind  maxima  and  their  signi¬ 
ficance  for  the  growth  of  nocturnal  inversions.  Bull.  Amer. 
Meteor.  Soc. ,  38^  283-290. 

Bleeker,  W. ,  and  J.  Andre,  1951:  On  the  diurnal  variation  of  pre¬ 
cipitation,  particularly  over  central  U.  S.  A.,  and  its  relation 
to  large-scale  orographic  circulation  systems.  Quart.  J.  Roy. 
Meteor.  Soc. ,  77, 260-271 . 

Bonner,  W.  D.  ,  1968:  Climatology  of  the  low  level  jet.  Mon.  Weather 
Rev.  ,  96,  833-850. 

_ ,  and  Jan  Paegle,  1970:  Duirnal  variations  in  boundary 

layer  winds  over  the  south-central  United  States  in  summer. 

Mon.  Weather  Rev.,  98 ,  735-744. 

,  S.  Esbensen  and  R.  Greenberg,  1968:  Kinematics  of  the  low- 
level  jet.  J.  Appl .  Meteor.,  7,  339-347. 

Buajitti,  K.  and  A.  K.  Blackadar,  1957:  Theoretical  studies  of  duirnal 
wind- structure  variations  in  the  planetary  boundary  layer.  Quart. 
J.  Roy.  Meteor.  Soc.,  8^,  486-500. 

Carl,  D.  M. ,  T.  C.  Tarbell  and  H.  A.  Panofsky,  1973:  Profiles  of  wind 
and  temperature  from  towers  over  homogeneous  terrain.  J.  Atmos. 
Sci. ,  30,  788-794. 

Curtis,  R.  C.  and  H.  A.  Panofsky,  1958:  The  relation  between  large- 
scale  vertical  motion  and  weather  in  summer.  Bull.  Amer.  Meteor. 
Soc.  ,  39.,  521-531  . 

Deardorff,  J.,  1974:  Three-dimensional  numerical  study  of  the  height 
and  mean  structure  of  a  heated  planetary  boundary  layer.  Bound. 
Layer  Meteor. ,  7,  81-106. 


j 


163 


Djolov,  G.  D. *  1973:  Modeling  of  interdependent  duirnal  variation  of 
meteorological  elements  in  the  boundary  layer.  Ph.D.  Thesis, 
University  of  Waterloo,  Waterloo,  Ontario. 

Dresden,  A.,  1930:  Solid  Analytical  Geometry  and  Determinants. 

Chapman  and  Hall  ,_Limi ted,  London,  118-121. 

Dutton,  J.  A.,  1976:  The  Ceaseless  Wind.  McGraw-Hill,  256-274. 

Feussner,  K.  and  P.  Dubois,  1930:  Trubungsfaktor ,  precipitable  water, 
staub.  Gerlands  Beitr.  Geophys.,  27,  132-175. 

Gerrity,  J.  D. ,  1976:  A  physical  numerical  model  for  the  prediction 
of  synoptic-scale  low  cloudiness.  Mon.  Weather  Rev.,  95,  261-282. 

Hewson,  E.  W. ,  1943:  The  application  of  wet-bulb  potential  temperature 
to  air  mass  analysis.  Quart.  J.  Roy.  Meteor.  Soc.,  63 ,  323-335. 

Holton,  J.  R. ,  1967:  The  diurnal  boundary  layer  wind  oscillation  above 
sloping  terrain.  Tellus,  19,  199-205. 

Jenne,  R.  L.  and  W.  Spangler,  1976:  USAF  average  elevation  data  -  one 
degree  and  5  minute  (Boulder,  Colorado:  National  Center  for  Atmos¬ 
pheric  Research). 

Kincer,  J.  B.,  1916:  Daytime  and  nighttime  precipitation  and  their 
economic  significance.  Mon.  Weather  Rev.,  44,  628-633. 

Lindzen,  R.  S.  and  H.  L.  Kuo,  1969:  A  reliable  method  for  the  numerical 
integration  of  a  large  class  of  ordinary  and  partial  differential 
equations.  Mon.  Weather  Rev.,  97,  732-734. 

Liou,  K.  N.,  1979:  Dept,  of  Meteor.,  Univ.  of  Utah  (personal  communi¬ 
cation)  . 

Lumley,  J.  and  H.  Panofsky,  1964:  The  Structure  of  Atmospheric  Turbu¬ 
lence.  John  Wiley  and  Sons,  Inc.,  New  York. 

McClatchey,  R.  A.,  R.  W.  Fenn,  J.  E.  Selby,  F.  E.  Volz  and  J.  S. 

Garing,  1971  :  Optical  properties  of  the  atmosphere.  3rd  ed. 
AFCRL-72-0497 . 

McDonald,  J.  E.,  1960:  Direct  absorption  of  solar  radiation  by  atmos¬ 
pheric  water  vapor.  J.  Meteor,  17,  319-328. 

Means,  L.  1.,  1952:  On  thunderstorm  forecasting  in  the  central  United 
States.  Mon.  Weather  Rev.,  80,  165-189. 

_ ,  1944:  The  nocturnal  maximum  occurrence  of  thunderstorms  in 

the  midv'tjstern  states.  Dept,  of  Meteor.,  Univ.  of  Chicago,  Misc. 
Rep.  No.  16,  37  pp. 


u. 


,  k - i.^‘IL 


1. 


m 


164 


Monin,  A.  S.  and  A.  Yaglom,  1971:  Statistical  Fluid  Mechanics; 

Mechanics  of  Turbulence.  MIT  Press,  Cambridge,  Mass.,  373-383. 

Ooyama,  Katsuyuki,  1957:  A  study  of  diurnal  variations  of  wind  caused 
by  periodic  variation  of  eddy  viscosity.  Final  Report,  AF  19(604) 
-1368,  College  of  Engineering,  New  York  University,  80-135. 

Orlanski,  I.,  B.  B.  Ross  and  L.  J.  Pol  insky,  1974:  Diurnal  variation 
of  the  planetary  boundary  layer  in  a  mesoscale  model.  J.  Atmos. 

Sci . ,  27,  965-989. 

Paegle,  J.,  1978:  A  linearized  analysis  of  diurnal  boundary  layer  con¬ 
vergence  over  the  topography  of  the  United  States.  Mon.  Weather 
Rev.  ,  106,  492-502. 

_ ,  and  G.  E.  Rasch,  1973:  Three-dimensional  characteristics 

"  of  diurnal ly  varying  boundary  layer  flows.  Mon.  Weather  Rev., 

101,  746-756. 

_ ,  and  D.  W.  McLawhorn,  1973:  Correlation  of  nocturnal 

thunderstorms  and  boundary  layer  convergence.  Mon.  Weather  Rev., 
101,  877-883. 

_ ,  W.  G.  Zdunkowski  and  R.  M.  Welch,  1976:  Implicit  differ- 

encing  of  predictive  equations  of  the  boundary  layer.  Mon .  Weather 
Rev. ,  104,  1321-1324. 

Paltridge,  G.  W.  and  C.  M.  R.  Platt,  1976:  Radiative  Processes  in 

Meteorology  and  Climatology.  Elsevier  Scientific  Publishing  Co., 
New  York,  60-66. 

Perkey,  D.  J.  1977:  The  effect  of  geostrophic  vs.  non-geostrophic 

initial  winds  on  quantitative  precipitation  forecasts.  Presented 
at  the  Third  Conference  on  Numerical  Weather  Prediction  of  the 
A.  M.  S.,  Omaha,  Nebraska. 

Phillips,  H.,  1962:  Zur  theorie  des  tagesganges  der  temperatur  in  der 
bodennaken  atmosphare  und  in  ihrer  unterlage.  Z.  Meteorol . ,  16,  5. 

Pielke,  R.  A.  and  Y.  Mahrer,  1977:  Mesoscale  model  predictions  forced 
by  geographic  features.  Presented  at  the  Third  Conference  on 
Numerical  Weather  Prediction  of  the  A.  M.  S.,  Omaha,  Nebraska. 

Pitchford,  K.  L.  and  J.  London,  1962:  The  low  level  jet  as  related  to 
nocturnal  thunderstorms  over  midwest  United  States.  J .  Appl . 
Meteor. ,  ]_,  43-47 . 

Plate,  E.  J.  1971  :  Aerodynamic  Characteristics  of  Atmospheric  Boundary 
Layers .  USAEC  Division  of  Technical  Information  Extension,  Oak 
Ridge,  Tennessee. 

Roewe,  D.  and  K.  N.  Liou,  1978:  Influence  of  cirrus  clouds  on  the  in¬ 
frared  cooling  rate  in  the  troposphere  and  lower  stratosphere. 

J.  Appl.  Meteorol.,  17,  92-106. 


I 


Rogers,  R.  R.,  1976:  A  Short  Course  In  Cloud  Physics.  Pergamon  Press 
Inc. ,  New  York.  ”  ' 

Shuman,  F.  G.,  1957:  Numerical  methods  in  weather  prediction:  II. 
smoothing  and  filtering.  Mon.  Weather  Rev.,  85,  357-361. 

Snyder,  V.,  1914:  Analytic  Geometry  of  Space.  H.  Holt  and  Company, 

New  York. 

Spenser,  J.  W. ,  1971:  Fourier  series  representation  of  the  position  of 
the  sun.  Search,  2,  172. 

Stanley,  D.  0.  and  G.  M.  Jurica,  1970:  Flux  emissivity  tables  for 

water  vapor,  carbon  dioxide  and  ozone.  J.  Appl .  Meteor.,  9_,  365- 
372. 

Wallace,  J.  M.,  1975:  Diurnal  variations  in  precipitation  and  thunder¬ 
storm  frequency  over  the  conterminous  United  States.  Mon.  Weather 
Rev. ,  103,  406-419. 

Webb,  E.  K. ,  1970:  Profile  relationships:  the  log-linear  range  and 
extension  to  strong  stability.  Quart.  J.  Roy.  Meteor.  Soc.,  96, 
67-90. 

Welch,  R.  M.,  J.  Paegle  and  W.  G.  Zdunkowski,  1978:  Two-dimensional 
numerical  simulation  of  the  effect  of  air  polution  upon  the  urban- 
rural  complex.  Tellus ,  30,  136-150. 

Wexler,  H.  1961:  A  boundary  layer  interpretation  of  the  low  level  jet. 
Tellus,  13,  368-378. 

Yamada,  T.  and  G.  Mel  lor,  1975:  A  simulation  of  the  wangara  atmos¬ 
pheric  boundary  layer  data.  J.  Atmos.  Sci.,  32,  2309-2329. 

Yu,  T.  W.,  1977:  Notes  on  a  turbulence  parameterization  of  the  atmos¬ 
pheric  boundary  layer.  J.  of  the  Meteor.  Soc.  of  Japan,  55,  617- 
622. 

Zdunkowski,  W.  G.  R.  M.  Welch  and  J.  Paegle,  1976:  One-dimensional 
numerical  simulation  of  the  effects  of  air  polution  on  the  plane¬ 
tary  boundary  layer.  J.  Atmos.  Sci.,  33,  2399-2414. 

Zdunkowski,  W.  G. ,  J.  Paegle  and  J.  P.  Reilly,  1975:  The  effect  of 
soil  moisture  upon  the  atmospheric  and  soil  temperature  near  the 
air-soil  interface.  Arch.  Met.  Geoph.  Biokl.,  24,  245-268. 


NAME 


David  Wilton  McLawhorn 


BIRTHPLACE 
BIRTHDATE 
HIGH  SCHOOL 

UNIVERSITIES 


DEGREES 


PROFESSIONAL  AND  HONORARY 
ORGANIZATIONS 


Ayden,  North  Carolina 

24  February,  1943 

Winterville  High  School 
Winterville,  North  Carolina 

East  Carolina  University 
Greenville,  North  Carolina 
1961-1965 

The  Pennsylvania  State  University 
University  Park,  Pennsylvania 
1965-1966 

University  of  Utah 
Salt  Lake  City,  Utah 
1971-1972 
1976-1979 

Bachelor  of  Arts  in  Mathematics,  1965 
East  Carolina  University 

Bachelor  of  Science  in  Meteorology, 
1966 

The  Pennsylvania  State  University 

Master  of  Science  in  Meteorology, 

1972 

University  of  Utah 

Sigma  Xi 
Pi  Mu  Epsilon 
Chi  Epsilon  Pi 
Phi  Kappa  Phi 

"The  effects  of  Weather  on  propaga¬ 
tion  of  electromagnetic  energy,"  Air 
Defense  Command  Communication  and 
Electronics  Digest,  15-70,  1970. 


PUBLICATIONS 


"Nocturnal  Thunderstorms  and  Bound¬ 
ary  Layer  Convergence",  (Masters 
thesis.  University  of  Utah,  1972), 

147  pages. 

"Nocturnal  Thunderstorms  and  Bound¬ 
ary  Layer  Convergence",  1973,  with 
Jan  Paegle.  Three  Dimensional  Char¬ 
acteristics  Of  Diurnal ly  Varying 
Boundary  Layer  Flows  (Technical  re¬ 
port  to  the  National  Science  Founda¬ 
tion),  129-275. 

"Correlation  of  Nocturnal  Thunder¬ 
storms  and  Boundary-Layer  Conver¬ 
gence,"  1973,  with  Jan  Paegle. 

Monthly  Weather  Review,  101,  877-883. 


